a = a <= = 3 - 
oval BU SSSRERSESRSEHSS SS SeAcABSSSEeeeS Co eae eee ee Pe cc lc cChLcrLhCLULC Cc crc Cc CC 


& CSREES RES 
MSI 204 1s ote 
ERR U Sh = CORE RDASOS ESS LRRERE SERS SESSA BEES ESS SES DDOREEEe 
Beha EE —oo eee eee 


Be sk S SERRE AREHS CHER CKRETE SR ROR ASS ttt 
ltt phe ()nen t niversityv. Ht BR RRSSBRRERSS OT ee eee eae 
se FRRFRSRRER SRS RREEE aeccegea an sana 
| EGER STR GAMES OUIENICES /Technoto; +B SSBSRERBESS955o 2080 Urea Et 
ie ie aie ek id i ee ee ede 
iL Li4 cere eee eee eee 
tT LMS D004 Mathematical Models ana (Methods) titi t it i  P to 
Prtlit iii titi iti Peet LL LL 


Rh Guaene? 


* : 


. 
ti 
SoC ewes 


T7777 s/h hhh 
RBESEEGR’ o. MESSE SNase 


oo 6o8e 
a 


wey ERY eh rs mes ~% - re 
ee w Le 
se hn colt wa el a 


oo 
o 
ci 
Saeen. “SR 


_ 


- 
- a 
ar S 
Co 4 
Fitts LLL LA lh LhLCULDLcL Llc cc LCL CL 
SS ts, 4 


The Open University 


Mathematics/Science/Technology 
An Inter-faculty Second Level Course 
MST204 Mathematical Models and Methods 


Unit 18 
Polynomial approximations 


Prepared for the Course Team 
by Mick Bromilow 


The Open University 


The Open University, Walton Hall, Milton Keynes. 
First published 1982. Reprinted 1983, 1985, 1989, 1992, 1994, 1996. 
Copyright © 1982 The Open University. 


All rights reserved. No part of this work may be reproduced in any form, by 
mimeograph or any other means, without permission in writing from the publishers. 


Designed by the Graphic Design Group of the Open University. 
ISBN 0 335 14047 5 


Printed and bound in the United Kingdom by Staples Printers Rochester Limited, 
Neptune Close, Medway City Estate, Frindsbury, Rochester, Kent ME2 4LT. 


This text forms part of the correspondence element of an Open University Second 
Level Course. 


For general availability of supporting material referred to in this text, please write to: 
Open University Educational Enterprises Limited, 12 Cofferidge Close, Stony Stratford, 
Milton Keynes, MK11 1BY, Great Britain. 


Further information on Open University courses may be obtained from 
The Admissions Office, The Open University, P.O. Box 48, Milton Keynes, MK7 6AB. 


1.7 


MST204 Unit 18 


Contents 


Introduction 


5 


Study guide 


Basic theory 


1.1 Polynomials 

1.2 Taylor polynomials 

1.3 Error bounds for Taylor polynomials 
Summary of Section 1 


Two applications of Taylor polynomials 


2.1 Roots of equations: the Newton-Raphson method 
2.2 Taylor series methods for differential equations (Tape Subsection) 
Summary of Section 2 


The catenary (Television Section) 


3.1 Hyperbolic functions 
3.2 The rope problem 
Summary of Section 3 


Interpolation and integration 


41 Interpolation 
4.2 Methods of numerical integration 
Summary of Section 4 


End of unit exercises 


Appendix: Solutions to the exercises 


4 MST204 18.1 


Introduction 


Many methods of applied mathematics are based on the use of polynomial 
approximations, and in this unit we shall consider how some of these methods are 
derived and used to solve practical problems. In the next unit we shall make use of 
many of these results when solving differential equations numerically. 


We consider two ways in which we might use a polynomial to approximate a 
function f(x). If we know the value of the function and its first n derivatives at 
some point x = a, say, as in Figure 1, then we can determine the Taylor In this unit the word point 


polynomial of degree n about « which approximates f(x) near the point x = «. means a real number in the 
domain of the function. 
In many problems the values of the derivatives are not known but the values of f 


are given at n+ 1 points Xo, x1, X2,..., X, aS in Figure 2. In this case we can 
determine that polynomial of degree n which takes these values of f at 
X05X1,---,X,. Such a polynomial is called an interpolating polynomial. 


x 


y Xx 


Xo. xX) X30 x‘, x 


Figure 1. For Taylor polynomials 
the data is given at x = a. Figure 2. For interpolating polynomials the data is given at n + 1 points. 


Taylor polynomials form the basis for two important methods in this unit: a 
method for finding roots of equations and a method for finding the approximate 
solutions of differential equations. Further, in the television programme we look at 
how Taylor polynomials can be used to derive approximate’ relationships between 
variables where it is impossible to find an exact relationship. 


In the final section of the unit we use interpolating polynomials to derive methods 
for determining the approximate values of integrals, such as the trapezoidal 
method and Simpson’s method. 


Study guide 


Section 1 of this unit revises the determination of Taylor polynomials and is 
central to the rest of the unit. Subsection 2.1 should be studied before Section 3, 
while Section 4 can be studied at any time after Section 1. Subsection 2.2, the tape 
subsection, can be studied at any time after Section 1, and is an important 
introduction to the next unit. Section 5 contains some extra exercises which you 
may attempt if you need more practice on the techniques introduced in this unit. 


Before watching the television programme you should have studied Subsection 2.1, 
and Section 3 up to Subsection 3.3. The television programme concerns the study 
of freely hanging ropes and cables. 


| Basic theory 


~ 1.1. Polynomials 
A polynomial is a function, p(x) say, which can be expressed in the form 
P(X) = ag + ayx + anx* + °° +,x" 


where n is a non-negative integer and do, a;, ... , da, are real numbers called the 
coefficients of the polynomial. 
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If n is the highest power of x for which the coefficient a, is non-zero, p(x) 1 is said to 
be a polynomial of degree n. Hence the polynomial 


it Jen 


is a polynomial of degree 3. 


Polynomials of low degree have acquired special names. 


quintic 


Thus the polynomial 1 +7x + 4x° is referred to as a cubic polynomial or simply a 
cubic. 


The advantage of using polynomials is that they are very easy to manipulate 
algebraically. For example, if p(x) is given by 


p(x) = 1+3x + 4x? + x? 

and q(x) is another polynomial given by 
q(x) = 3 + 2x — x?, 

then we can, for example, add the two polynomials to get another polynomial 
p(x) + q(x) = 4 + 5x + 3x? + x? 

Polynomials are also very easy to differentiate or integrate. To differentiate 
p(x) = 1+ 3x + 4x? + x?, 


we simply differentiate each term separately to obtain 
# ix) = 3 + 8x + 3x’, 

while integrating p(x) gives 
ove) ax = fa + 3x + 4x? + x?) dx 


3 4 1 
=C+xta toe toe 


where C is an arbitrary constant. 


Example 1 


The function tan x may be approximated near x = 0 by the polynomial 


Use this approximation to estimate 


6.5 
{ tan x dx. 


0 


Solution 


0.5 0.5 2 
| tanxdx ~ ie be + 4+ — a) a 
0 15" 


2 
=15% + 45° |, 


= 0.130556 to six decimal places. 


These definitions are also 
given in Subsection 3.1 of the 
Handbook. 


p(x) is the 5th Taylor 
polynomial about 0 for tan x. 
Taylor polynomials will be 
discussed in Subsection 1.2. 
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(The true solution, using the integration tables in Subsection 7.1 of the Handbook, 


gives 0.5 O.3 
{ tan x dx = —log, cos :| 
8) 0 


= 0.130584 to six decimal places. 
Hence our approximation method gives a result which is correct to four decimal 
places.) 


When evaluating a polynomial for a particular value of x, the most efficient 
method to use is called nested multiplication. This method will also usually reduce 
rounding errors in the calculations. To see how the method works we consider the 
polynomial 


p(x) = 1 + 3x + 7x? + 4x°. 


If, for convenience, we reverse the order of the terms in the polynomial so that we 
write 


p(x) = 4x? + 7x? + 3x +1, 

then the nested multiplication method can be used to evaluate p(x) as 
p(x) = ((4x + 7)x + 3)x + 1. 

The steps in the calculation are as follows. 


(i) Enter 4 into the calculator. 


(ii) Compute 4x43 On my calculator the 


sequence of operations, given 


ili Compute (4x + 7) x x 4+ 3. that x is stored in the memory 
( Se and can be obtained using the 
(iv) Compute ((4x + 7)x +3)x x41. RCE keys 
This sequence of operations can be expressed using a recurrence relation as A be 4] [x] [RCL] Fe [+] +] 7] a 

Up = 4 


Uy = XUy + 7 


‘lg 
ellie’ 
ad Oa 
steal 
EE 


[=| 
=] 


Dx] Ex] 4 


uy = Xu, + 3 
uz = Xu, + l, 


and u3 gives the value of the polynomial at x. 


Example 2 
Evaluate p(x) = 4x° + 7x? + 3x + 1 at x = 1.27. 


Solution 
Using the recurrence relation, we have 
Up = 4 
uy = 1.27ug + 7 = 12.08 
u, = 1.27u, + 3 = 18.3416 
uz = 1.27u, + 1 = 24.293 832. 


Notice that this method of evaluating the cubic requires only 3 multiplications and 
3 additions. 


The following procedure box gives the recurrence relation for evaluating a 
polynomial of degree n. 
Nested multiplication 


To evaluate 


P(X) = A,X" + Aq—1X""* + +++ + a,x + do 


at a particular value of x, the nested multiplication method uses 
the recurrence relation 


Up4 1 = XU, + An—y-1 


where up = a, 
and u, gives the value of p(x). 
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Exercise 1 
This exercise concerns the polynomial 


Bale ae +38 4 x? — Tx + 2 


(i) What is the derivative of p? 

(ii) What is the integral of p? 

(iii) Evaluate p(x) for x = 0.29 using nested multiplication. 
(iv) Evaluate the derivative of p at x = 0.77. 


(Solution on p. 50] 


1.2 Taylor polynomials 


There are many functions such as sin x, e*, log, x and so on which we have to 
evaluate for a given value of x when carrying out computer or calculator 
computations. The traditional method, which you may have used, was to look up 
the required value of sin x in a table of sines. It is however not feasible for 
computers or calculators to store tables of numbers, since it would use up a 
relatively large amount of computer time to locate the required value. What 
happens in practice is that we approximate functions like the sine function by 
polynomials, so that for a given value of x we approximate sin x by the value of 
the polynomial. | 


There are several types of polynomial approximation, but traditionally the most 
popular are the Taylor polynomials. Suppose that we want to determine a_ 
polynomial p(x) of degree n which approximates some given function f(x) near 
x = 0. We write down p(x) as 

P(X) = do + aX + anx* +--+ + a,x". 


To derive the Taylor polynomial, we require p to be such that 


p(0) =f (0), 
p'(0) =f'(0) 


p’(0) =f" (0), 


p (0) = f(0). 
This uniquely determines p, since 
p(0)=ao, _—so that ay = f (0); 
p'(0)=a,, so that a, = f'(0); 


(0) =k!a,, so that a, _ = 
Hence 
2 n 
p(x) =f) + xf'O) + FF'O) + += FO) (1) 


This polynomial will normally give very good approximations to f near x = 0. 
(However, there are some functions, such as |x|, which cannot be approximated in 
this way since their derivatives do not exist at x = 0.) The above polynomial is 
known as the nth Taylor polynomial about 0 for the function f. 


Here are two examples of Taylor polynomials. 


Example 3 
Determine the 3rd Taylor polynomial about 0 for the function 


faye. 


This derivation is also given in 
the Taylor Series units of 
M101 and M8283. 


You have already used one 
Taylor polynomial in Section 
3 of Unit 7. See also Section 4 
of the Preparatory Booklet, in 
which the alternative name 
nth Taylor approximation is 
used. 


8 MST204 18.1 


Solution 

Differentiating this function is easy, since 
d 
—e=e, 
dx 

so that 


fW=e, f"(x)=e and f"(x)=eN%. 
At x = 0, e& =1,s0 that 


IO =F OF =F DO) Hf = t 
Hence the 3rd Taylor polynomial about 0 approximating f(x) = e* is obtained by 
substituting into Equation (1), giving 
eee 
2 
In Figure 1 we have plotted the graphs of e* and p(x). Note that the two graphs 
are indistinguishable near x = 0. 


p(x) =1+x+ 


y=ltxthx2+ tx) 


Figure 1. Graphs of e* and p(x) =1+x+ =a + ax 


Example 4 
Consider the function 


f (x) = (1 + xy’, 


where n is a positive integer. This is a polynomial of degree n and we would 
normally use the Binomial Theorem to obtain the expansion. However, we shall The Binomial Theorem is 


; ; | ; given in M101 and MS283, 
obtain the same polynomial expansion by looking at Taylor polynomials. Set actettian tatthe 


The derivatives of f(x) are Preparatory Booklet. 


f(x) =n(l + xy"? 
f"(x) = n(n — 1)(1 + xy"? 


fO(x) =n! 
All higher derivatives are zero. 


At x = 0, (1 + x) = 1 for any value of k, so that 
fO=1, f\O=n, f"O)=n(n—-1), ..., fMO)=n! 
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Hence, using Equation (1), we have the nth Taylor polynomial about 0 for 
f(x) = (1+ x)" as | 


— 1 —1)(n-—2 
p(x) =1+nx + ee + men ens ee. 
Here p(x) is exact because (1 + x)" is a polynomial of degree n, and we can write 


Se yO A 


n 
1 ee | 
(i + x) 5 eae 31 
In an exercise at the end of this subsection you will be asked to determine a 
Taylor polynomial approximation for (1 + x)!/. 


Very frequently we need approximations near some point x = a where a is not 
necessarily equal to zero. Suppose that we look for a polynomial of the form 


P(X) = apo + ay (x — a) + an(x — a)? +--+ +.,(x = a)" 
such that 


p(a) = f (a) 
p'(a) =f '(«) 


pa) = f(a). 


Such a polynomial is called the nth Taylor polynomial about « for the function f. 
Again it is not difficult to deduce that the coefficients of this polynomial are 


_f%@) 
a = k! 9 


Thus we may define these Taylor polynomials as follows. 


The nth Taylor polynomial about « approximating the function 
f (x) is defined by 


— 2 — 
plx) = f(a) + ( — a) f(a) + Spay + ++ 4 FP poyay 


To evaluate polynomials of the form 


P(x) = do + ay(x — a) + aa(x — a)? +--- +.a,(x — a)", 
the nested multiplication procedure can be adapted as 
Ur+1 = (x a ou, + An—p—1 


where uo = a, and u, gives the value of p(x). Thus polynomials in this form are as When evaluating them, do not 


easy to evaluate as are ordinary polynomials. multiply out (x — «)’ as 
x? — 2ax + a?, etc. 


Example 5 
Determine the 3rd Taylor polynomial about 1 approximating the function 
fRl=e. 
Evaluate this polynomial at x = 1.55. 
Solution 
For n = 0, 1, 2,... we have f(x) = e* and f™(1) = e. Thus, 
p(x) =e + (x — lhe + 4(x — 1)?e + 3(x — 1)%e 
=e{1+ (x — 1) + Hx — 1)? + &x — 1)*}. 


This approximation is compared with e* in Figure 2 (overleaf). In this case the 
graph of the polynomial is indistinguishable from that of e* near x = 1. 
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y y=e* 
10 
8 
6 
4 


yoe(ht («= 1) tog eo 


2 
—2 ~1 0 1 eS 


Figure 2. Graphs of p(x) = e{1 + (x — 1) + (x — 1)? + #(x — 1)°} and &. 


To evaluate p(x) at x = 1.55, compute (x — 1) = 0.55. Using nested On my calculator I stored 
multiplication to evaluate 1 + (x — 1) +4(x — 1)? + d(x + 1)3, we have 0.55 in the memory using the 
RCL key to recall it. Here is 
Uo = 0.166 666 67 © my key sequence 


uy = 0.55u9 + 0.5 = 0.591 666 67 


G 
ur = 0.55u; + 1 = 13254167 =] [x] [Rex] FF] 11] SI 
us = 0.55u2 + 1 = 1.7289792. RCL =] 


Hence p(1.55) = eu3 = 4.699 85 to five decimal places. This is not a bad 
approximation to the correct answer, which is 4.71147 to five decimal places. 


Exercise 2 


Determine the 2nd Taylor polynomial about 0 approximating f(x) = (1 + x)'/2. Evaluate 
this polynomial at x = 0.44 and compare your answer with the true solution. 


Exercise 3 
Show that the nth Taylor polynomial about 0 approximating f (x) = log, (1 + x) is 
a | n+1 x" 
eee 
n 


p(x) = x — 4x? + 43 - 


[Solutions to Exercises 2 and 3 on p. 50]. 


1.3 Error bounds for Taylor polynomials 


Wherever we use these approximations it is wise to obtain some estimate of the 
error in the approximation at a given value of x. That is, we would like to know 
something about the magnitude of the error function ¢(x) defined by 


e(x) = p(x) — f(x). | (2) 


For example, the error function for the third Taylor polynomial about 0 
approximating e*, given in Example 3, is 


e(x) =14+x+4x7 +309 -—e&. 
(x) has been plotted in Figure 3. 
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\ 
\ 
| 
| 
y = -0.187x4 | 


Figure 3. Graph of e(x) = 1 + x + 4x? + 4x3 — &, 


For comparison, the graphs of the functions +0.187 x* are also shown. (The 
reason for choosing these particular functions will become clear when you read 
Example 6 on pp. 12-13.) 


The behaviour of e(x) near x = 0 in the above example should be expected because 
we have required that the first three derivatives of p(x) and e* should agree at 

x = 0. Hence «(0) = e’(0) = ”(0) = e”’(0) = 0, and so the first three Taylor 
polynomial approximations to ¢(x) are all zero. 


Suppose now that we have obtained the nth Taylor polynomial p(x) about « 
approximating some given function f(x), and we wish to use this polynomial p to 
approximate f at a particular value of x. We need some idea of what the 
magnitude of the error is at this value of x. The theorem that follows gives a 
method of doing this. It is usually stated not in terms of the error function &(x) 
but in terms of a function called the remainder, defined as the function r(x) in the 
formula | 


I(x) = p(x) + r(x). 


Comparing with (2), we see that the remainder and the error function are related 
by 


r(x) = —e(x). (3) 


Taylor’s Theorem 

Let f and its first n + 1 derivatives exist and be continuous at all 
points between and including « and x. Then f(x) can be expressed 
as 


om 
f(x) =f (a) + (x — a) f(a) +o pa+ <e 


oe * ea pa) 4 r(x) 


where r(x), the remainder, satisfies 
= (x ees ar bai 
r(x) = ep senee (cx) 


for some number c,. between « and x. 


The expression (4) is often called the nth Taylor expansion of f (x) about a with 
remainder. The infinite sum 


£0) =F (2) + (= a)f'(a) +S pray ¢ FIM yma g 5 


The requirement that the 

n+ 1 derivatives of f(x) 
should exist excludes functions 
such as |x| or x!/3 over an 
interval which includes 0. 


I have used c, here rather 
than c to indicate that c, 
depends on the value of x. 
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is called the Taylor expansion (or Taylor series) of f(x) about «. Its significance is 
that, if we know for some particular value of x that the expression for the 
remainder, 


(x — a)" * 1 
(n + 1)! 


approaches zero as n becomes large, then we can conclude that (for that value of 
x) the sum of the terms in the Taylor series converges to f(x). Although we will 
not be going into the question of convergence, the use of the Taylor series 
expression for a function is so widespread in applied mathematics that it is 
important to remember expression (5). 


x) ee ee) 


Rather than prove Taylor’s Theorem I will just show what happens in the 
particular case when n = 0. In this case Taylor’s Theorem states that if f and f’ 
exist and are continuous for all points between and including « and x, then 


f(x) =f (a) + (& — a)f (cx) | (6) 
where c,, lies between « and x. 
Rearranging Equation (6) gives 


I(x) ~f(%) a =f) 


— 


F'(ex) = 


Graphically this can be interpreted as saying that between « and x there must be a 
point c, for which the slope of the graph (f‘(c,)) is the same as the slope of the 
(dashed) line joining the points («, f(«)) and (x, f(x)) in Figure 4. 


(x, f(x)) 


se | 
Sau slope of this 
line is 


I(x) FQ) 


See 


Figure 4. Graph illustrating Taylor’s Theorem for n = 0. 


This special case with n = 0 has an alternative name. It is often called the Mean 
Value Theorem for derivatives. 


Since the remainder and the error are equal and opposite, we can use Taylor’s 
Theorem to estimate the error function eé(x). The following example shows how 
this can be done. 


Example 6 
(i) Verify that for all c satisfying 0 < 


C 
(ii) Deduce that for all x satisfying 0 < 
polynomial about 0 for e* satisfies 


< 1.5 we have |e’| < e' > = 4.4817. 
x < 1.5 the error of the 3rd Taylor” 


1 
|e(x)| < Ax*, | where A = —e!°>. 
24 
Solution 
(i) Since e* is always positive we have |e‘| = e, and since c < 1.5 and the 
function is increasing we have e° < e!:> = 4.4817. 


It follows that le‘| < 4.4817... as required. 
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(11) Taylor’s Theorem, with f(x) = e*, gives 
e*=1+x + 3x? + dx? + r(x), where 
r(x) = ax* fO(c,) 
= ax" exp (c,). (7) 


Since 0 < x < 1.5, it follows that 0 < c, < 1.5 and hence, by part (i), that 
lexp(c,)| < e*°. 


Using (7) we deduce that 


Ir(x)| < zax* el 


= Ax*, where A = 4e!-5 ~ 0.187. 
Since ¢(x) = —r(x) (Equation (3)), we have |e(x)| = |r(x)|, and hence 
|e(x)| < Ax* 


as required. 
This result can also be written 
—Ax* < &(x) < Ax’. 


Figure 3 (page 11) shows that these inequalities are indeed satisfied for 
O <x < 1.5, as the theorem states. 


The procedure for determining the error bound used in Example 6 is formalized in 
the following box. 


Error bounds for Taylor polynomials 


1. Let p(x) be the nth Taylor polynomial about « approximating 
a given function f(x) which satisfies the conditions of Taylor’s 
Theorem. We want to find a bound on the error function 


e(x) = p(x) — f(x) 


for all x satisfying a < x < b, where a, b are given numbers 
such thata<aanda<b. 


2. Determine a value for M such that 
te te = M 


for all c satisfying a<c<b. 


3. Taylor's Theorem gives 


(x — a)"*? mare 
é(x) = ———___ C 
where Cy lies between « and x, and it follows since a < c, <b, 
that 


|x hss aoe 
(n+ 1)! 


This is the required error bound. 


le(x)| < 


Exercise 4 
(i) Verify that for 0 < c < 0.5 we have |1 + c| < 1.5. 
(ii) Compute an error bound for the 3rd Taylor polynomial about 0 approximating 


(1 + x)°, valid for all x satisfying 0 < x < 0.5. (Note that you are not asked to 
compute the 3rd Taylor polynomial itself.) 


Exercise 5 ; 
(i) Verify that |sin c| < 1 for all c. 
(ii) Determine the Sth Taylor polynomial about 0 approximating sin x. Compute an 


error bound for this approximation at x = 1, using Taylor’s Theorem. 


[Solutions to Exercises 4 and 5 on p. 50] 


Remember that exp (...) 
means e to the power of... . 
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Summary of Section 1 | 
1. A polynomial is a function, p(x) say, which can be expressed in the form 
p(x) = do + ax + anx* +°°* + a,x" 


where n is a non-negative integer and do, a;,..., a, are real numbers called the 
coefficients of the polynomial. If a, is non-zero, p is called a polynomial of degree n. 


2. To evaluate p(x), written as 
P(X) = qx" + Ani X"~* +°°+ + 4X + Ao, 


for a particular value:of x, the nested multiplication method uses the recurrence 
‘relation 


Urp+1 = XU, + An-r-1 


- where uy = a, 
and u, gives the value of p(x). 


If p(x) is of the form 
p(x) = a,(x — a)" + a, 1(x — a)" 4+ "++ + a(x — @) + ao 
then the recurrence relation is | 
Up +1 = (X — Xu, + An—y—1 


where, as before, up = a, and u, gives the value of p(x). 


3. The nth Taylor polynomial about « approximating fis given by 


@ — a), (= a) 
set Na) tok Sf) 


p(x) =f (a) + (x — a) f(a) + 


4. Taylor’s Theorem. Let f(x) and its first n + 1 derivatives exist and be 
continuous at all points between and including « and x. Then f(x) can be 
expressed as 


fe) =f@+e-—asia+® me Pia = 
(x —a)" Ce | aren 
a eee 


; x-2@ n+1 
where c, lies between « and x and the expression aan f"* (c,) is the 
remainder. (n+ 1)! 


Hence the error function ¢(x) can be written as 
e(x) = p(x) — f(x) 


= (x cm ay = 
= Bree. 1c). 


If a number M can be found such that 
If"* Y(c)}] <M 


for all c satisfying a < c < 5, then it follows that 


This is an upper bound on the error. 
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2 Two applications of Taylor polynomials 


2.1 Roots of equations: the Newton-Raphson method 


In the course of solving mathematical problems it is often necessary to find the 
values of x for which the equation 


f(x) =0 


is satisfied, where f is some known function. I shall refer to them as the roots of 
the equation f(x) = 0. 


One way of determining the roots of an equation is to sketch the graph of 
y =f (x) 

and to determine those values x for which the curve cuts the x-axis. 

For example, to find the roots of the equation 
f(x) = 2x3 — 5x? — 5x + 10 =0 

we could sketch the graph as in Figure 1. 


y= 22°-Sx*-5x +10 


Figure 1. Sketch of y = 2x> — 5x? — 5x + 10. 


15 


They can also be called the 
solutions of the equation. You 
have already met these words, 
used interchangeably, in 

Unit 5, Complex numbers. 


Curve sketching is reviewed in 
Subsection 5.3 of the 
Handbook. 


Looking at Figure 1 we can see that there are three values of x at which the graph 


cuts the x-axis. These are approximately — 1.5, 1.2 and 2.8. 


In very rare cases the graph may only touch the x-axis but not cross it, so that the 


graph may be as in Figure 2. 


Figure 2. Graph of a function 
which touches but does not 


It can be difficult to determine roots in this instance, since with the naked eye we cross the x-axis. 
cannot detect whether there is one root, two roots or no roots as in the following 


close-ups. 


(a) one root (b) two roots (c) no roots 
Figure 3. Possible close-up views of the graph in Figure 2. 
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However in this unit we will assume that the graph is of the type given in Figure 1. 


Another method of finding roots, which you have met before, is the formula — See Subsection 3.1 of the 
method for quadratic equations such as Handbook. 


x? — 13x +7=0. 
Here are some equations for which the roots are not so easy to compute. 
(i) x —2sinx =0 
(ii) x3 + 13.5x? + 40x + 16.67 =0 
(111) xe *—0.1 =0 
What we want is a method which could be used to solve any of these problems. 


The method we will use is an application of Taylor polynomials and is called the 
Newton-Raphson method. 


The Newton-Raphson method 


This method determines the roots of an equation using a recurrence relation. We The Newton-Raphson method 
begin by supposing that there is a root of the equation f(x) = 0 near some point was introduced in the Taylor 
Xo. Now the first Taylor polynomial about x, approximating f(x) is given by Series units of M101 and 


MS283. 
P(x) =f (Xo) + (* — Xo) f'(Xo). (1) 


This approximation is often called the tangent approximation to f(x) at xo because 
it is the equation of a straight line which is tangential to the curve y = f (x) at Xo, 
as in Figure 4. 


the tangent 
approximation 
at Xo, 


p(x) =f(X9) + (x — Xo) S'( Xo) 


Figure 4. The tangent approximation to f(x) at Xo. 


In Figure 4, xo is not a particularly good approximation to the root of f(x) = 0, 
but we can use the tangent approximation at Xo to improve this. The graph of the 
tangent approximation cuts the x-axis at the point x, such that 


p(x1) = 0. 


Since p is an approximation to f, the root, x,, of the equation p(x) = 0 is likely to 
be a close approximation to the root of f(x) = 0. 


To calculate x, we use Equation (1), which tells us that p(x,) = 0 when 
f (Xo) + (%1 — Xo) f'(Xo) = 0. 


Hence rearranging gives 


f (Xo) A slightly different derivation 
ae | 2 of Equation (2) will be given 
X14 Xo ; : ( ) : ie 

f'(Xo) in the television programme. 
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Example 1 
From Figure 1 we know that there is a root of the equation 


f (x) = 2x*® — 5x? — 5x + 10 =0 Remember to use nested 
multiplication to evaluate f(x). 
near X> = —1.5. 
Now f (—1.5) = —0.5 
and f(x) = 6x? — 10x — 5, giving 


f'(—1.5) = 23.5. , 
Using Equation (2) with x) = —1.5 gives 
f(-1.5) 0.5 
= -15- = —15 + — = —1.478 7234. 
“1 FEAL.) 7 ae 


We can see that x, is probably a much better approximation to the root than Xo, 
since f (— 1.478 723 4) = — 0.006 318 39 and this is considerably smaller in 
magnitude than f(—1.5) = —0.5. 


We need not stop here. We can compute the tangent approximation to fat x, to 
obtain an even better approximation, x,, to the root of f (x) = 0, as in Figure 5. 


the tangent 
approximation 
at Xo 


2) xy 


the tangent approximation 
at xX] 


Figure 5. The tangent approximation at x,. 


The equation for x, can be deduced from Equation (2) as 


f (x1) 
f'(x1) 
Using this value of x. we can compute x3 using another tangent approximation at 


X2, and so on. Each term in the sequence xo, X;, X3, ... is computed using the 
recurrence relation 


M2 = Xi = 


J (Xr) 


Xr+1 ar Xy f' ) 
r 


Example 2 


(3) 


In Example 1 we looked for a root of the equation 
2x> — 5x? — 5x +10=0 


near Xo = —1.5 and computed x, = — 1.478 723 4. Continuing with this example, 
we build up the following table using Equation (3). (You may get slightly different 
results with your calculator.) 
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: fA) 
f'@,) 


f (x,-) = 2x,? — 5x,” — 5x, + 10 f'(x,) = 6x,” — 10x, — 5 


Xr+i1 sa Xy 


0} —15 —0.5 23.9 — 1.478 723 4 
1 | —1.478 7234 — 0.006 318 39 22.906 971 — 1.478 4476 
2 | —1.478 4476 — 0.000 001 06 22.899 S19 — 1.478 4475 
3 | —1.478 4475 —2x 107° 


After only three iterations f (x3) is so small that we can conclude that 
x3 = —1.478 4475 1s a very good approximation to the root. — f (x3) depends quite strongly 


on which calculator is used 
and how the calculation is 


The precise value obtained for 


There are two possible criteria for stopping the sequence: 


(i) If f (x,) is very small, then x, is assumed to be very close to the root. nein plier ies aes d 
(11) If x,41 ~ x, then Equation (3) gives 
xX, ary ee ‘ X> 
X, 2X, — I(r) which is equivalent to J (%) ~ 0 


| f'() I )s = 
Hence provided that f’(x,) is not too large we can conclude that x, is again a good 
approximation to the root. 


If we apply these criteria to Example 2, we find respectively: 
(i) f (x3) = —2 x 10°°; 
(ii) X2 = —1.478 4476, x3 = —1.478 4475; 


and in either case we could conclude that there is no point in computing further 
terms in the sequence. 


We summarize the Newton-Raphson method in the following procedure box. 


The Newton-Raphson method 
Given an initial approximation Xp, close to the desired root of 


f(x) = 0. 


Generate a sequence Xo, X1, X2, ... of successively better 
approximations using the recurrence relation 


7) 
f'@,) 


This procedure can be terminated if 


Xr+1 =X, 


either: (i) f (x,) becomes sufficiently small 


or A Ss. 


To begin using the Newton-Raphson method we need an initial approximation Xp. 
If we do not have any information about x, the simplest thing to do is to sketch a 
graph of y = f(x), as we did for Example 1. The sketch need not be particularly 
accurate, and plotting half a dozen points will usually be sufficient to obtain crude 
approximations to the roots. In the following example, we suggest an alternative 
curve sketching method which sometimes can reduce the effort required to find a 
value for Xo still further. 


Example 3 
Find all the roots of 


x —2snx = 0 


using the Newton-Raphson method. Quote your answers to five significant figures. 
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Solution 


To use the Newton-Raphson method I need a value for x9. To obtain this I could 
sketch a graph of 


f(x) =x —2sinx 


and look for the roots. However, it is easier in this case to draw graphs of y = x 
and y = 2 sin x and to look for points at which the two graphs intersect, since at 
these points we have x = 2 sin x. 


We can deduce that there can be no roots if |x| > 2, since —1 < sin x < 1, so that 
—2 <2sin x < 2 and the only possible roots of x = 2sin x are in the interval 
—2<x <2. 


root root 


root 


Figure 6. Graphs.of y = x and y = 2 sin x for —2 < xe Zz 
From the graphs in-Figure 6 we can see that there are three roots. The first is at 
x = 0, while the other two roots are approximately —1.9 and 1.9. 
The Newton-Raphson recurrence relation is 
= fia) 
f'(x) 
In this case, f(x) = x — 2 sin x and f’(x) = 1 — 2cosx. Hence the recurrence 
relation is 


Xr+1 = Xp 


x, — 2sin x, 
Xp41 =X, —- ————. 
: " 1 —2cosx, 


With x, = 1.9 I obtained the following sequence, using my calculator. 
Xo = 19 
x, = 1.895 5059 
X27 = 1.895 4943 (f (x2) = —9 x 107°) 
x3 = 1.8954943 = (f(x3) = —4 x 107?°) 


After only three iterations the values of x, agree to eight figures and f(x3) is 
exceedingly small. 


Thus the root to five significant figures is 1.8955. 
Similarly the third root is at — 1.8955. Hence the three roots are 0 and +1.8955. 


In the above example the sequence converged very rapidly to the required root. 
The question we might ask is whether this will always be so. In the following piece 
of theory we will explain this rapid convergence and also gain some insight into 
how accurately we need to determine Xo, the initial guess at the root. 


Since we know that x = 0 is 
the exact value of a root, we 
do not need to use the 
Newton-Raphson method for 
this one. 


You may get slightly different 
results using your calculator. 


We know from basic 
trigonometry that 

sin (—x) = —sin x, so that if 
is a root of x = 2 sinx then so 
is —a. 
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The recurrence relation for the Newton-Raphson method given by Equation (3) 
can be rearranged as 3 


0 = f (x;) + (X,+ : Waa Xr)F '(X,). (4) 


Let us compare Equation (4) with the first Taylor expansion about x, with 
remainder (page 11). Assuming that f(x) satisfies the conditions of Taylor’s 
Theorem (with n = 1), we have 


f (x) = f (x,) 23 (x aS Xr) f (Xr) = (x poe Ky I Ae,) (5) 
where c,, lies between x, and x. 


If p is a root of the equation f(x) = 0, so that 


F (p ) = 0, 
then we can substitute p for x in Equation (5) to get 
0 =f (x,) + (0 — xf") + He — x) Fle) (6) 


where c, lies between x, and p. 
Subtracting Equation (6) from Equation (4) gives 
0 = Mee — pF’ O%) — He — x) F"CCp), 


X41 — P =p — noe = (7) 


Look at this last equation, (7), as this gives the secret of the success of the 
Newton-Raphson method. The left-hand side of the equation, x,+, — p, is just the 
error in the (r + 1)th approximation x,,, to the root p. Equation (7) relates this 
error to the error in x,. It shows that the error in x,,, is proportional to the 
square of the error in x,; soif the error in x, is fairly small, the error in x,4, will 
be much smaller. As soon as we are reasonably close to the root, the errors 
diminish very rapidly. 


In order to say more precisely when the method will converge, using Equation (7), 
we need to know some values for 


I"p) 
2f'(Xr) 


which are generally not available. However, we can say that provided the 
conditions of Taylor’s Theorem are met and the equation f’(x) = 0 does not have 
a root near p, then the sequence of iterates is guaranteed to converge for a suitable 
starting value x,. At the beginning of Subsection 2.1 I stated that we would 
assume that the function did not behave as in Figure 2. This figure shows a 
situation in which the equation f’(x) = 0 does have a root near p, so that the 
Newton-Raphson method would not work well. 


Exercise 1 
Find all three roots, to 5 significant figures, of the cubic equation 


x3 + 13.5x? + 40x + 16.67 =0 © 
using the Newton-Raphson method. Figure 7 shows a sketch of the graph. 
[Solution on p. 50] 


The results of this exercise will 
be used in Unit 24, Normal 
modes. 


If you are short of time, just 
calculate one of the roots. 
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yaxPHS 5x2 +400+16-67 7 
50 


40 


y(-l0)=-33.3 


Figure 7. Sketch of y = x? + 13.5x? + 40x + 16.67. 


Exercise 2 
(i) Show that the roots of the equation 
xe *—0.1=0 
are at the intersections of the graphs of y = e* and y = 10x, and hence find rough 
approximations to them. 
(11) Use the Newton-Raphson method to determine these roots to 6 decimal places. 
[Solution on p. 51] 3 


2.2 Taylor series methods for differential equations 
(Tape Subsection) 


The second application of Taylor polynomials that we look at is in the numerical 
solution of differential equations of the form 


y’ =m(x, y) with y(xo) given 
where m(x, y) is some given formula which may contain both x and y. 


For example, we may be asked to find an approximate solution to the differential 
equation 


ent y 
for values of x such that 0 < x < 1, given that y = 0 when x = 0. This is the 


example we shall use in the tape frames which follow. 


(In this particular case we could compute an analytical solution using the 
integrating factor method, but we use this simple example here to illustrate the 
numerical method.) 


You have already met one method, Euler’s method, in Unit 2 for solving such 
problems numerically. The methods we shall develop here are called Taylor series 
methods. In general Taylor series methods will give much more accurate results 
than Euler’s method. 


Start the tape now. 


You saw many differential 
equations of this form in 
Unit 2. 


The true solution is 
y=e-x-l. 
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Taylor polynomials 


Differential equation: y' = x+y (!) 


O 


Initial condition : y (0) 


PROBLEM : Use a 2nd Taylor polynomial approximation for 
y about O to get an approximate value for y ak x= 0'l 


p(x) = y(0) + acy'(0) + a oc* (0) (2) 


DATA : y (0) eae 
kg SSE 
oe 

~ - [or] 


Computing the second derivative 
The differential equation (1) is 
ye ery 


"= diu') = d (x+ 
Y Ay) a y) 


Louk ce SE 
Hence y'(0) = L__ 
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A new problem 
Differential equation : yy = (4) 


Initial condition : Y (0-1)= 0-005 


——= 
O.| 0.2 


PROBLEM: Use a 2nd Taylor polynomial about O to 
get an approximate value 4or y at 2 OD 


p(x) = y(ol)+ (x- 01) y’(01) + Ya(x-O')"y"(o4) (5) 


yo)? LJ 
oo Saas 
Wen: = 


Seuurion : p(oa) Saeemat 


4 | The general problem 
Differential equation : ry a (=, y) 


Initial condition : Uy (2<,) - x 


PROBLEM: Use a and Taylor Polynomial about =, to 
get an Approxmate value for 4 eG ac = ey) 


S G).> y (2¢,) + (2c - <-) uy! (2, ) + '/2 (nae jag =.) (6) 
where (sc, )= m (=<, 4 (*,)) (6a) 
snd uf (e,)= dm (=e (=) Ge) 


SOLUTION : p(=,,,)= 4(*,) + (20, 5 7.2¢,) ty 26) + Za. Feret” >,) y',)(7) 
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Simplification 
Write u, for Y(2c,), uy, for ui(ac,) and yu". for 4y"(2C,) 
Write h = 2¢,,,-2, as the step size 


Equation (7) becomes 


| P(-41) aie a hy, +/ah* ue 


where Yr =m (2, Yr) 


and yy = dev (2,.4,) 


Reworking the problem in 


Differential equalion : y mc + y 
Imitial condition : y (0) O 


Taylor series method of order 2 with 
~ f ! 2 wu 
p(x,) = y thy, + ah Yo 
where Uo mM (2c, ; U,) = 2c ,+ YU, = 
and u" dm (xc ee 
Mo am ( Yo) _ 


p(o'!) 


Notation for the approximate solution 
Given Y,, we want values for 4.,4,,44. °°, Yn 


Approximations : Given a zs Uy, we actually compute ee A Pata, & 
The Taylor series method of order 2 uses the recurrence relation 
Live - 7 +hy, = any, 

where Y- = m [ oer ved 


andY’ = dm (2,,Y,) 
osc 


a) 
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© Atypical problem 
Compute the approximate solution bo the diPPerential 
equation 


é 


ym ery, 


with (8) - O, at «sO 6-2-1] using the Taylor 
secies method of order 2. | 


er Ha Vet Rs + a Y 
0} 0 0 0 


0.005 
f/-0.1-1 “0:005 0.105 1.105 0.021025 
2} 0.2) 6021025 
31-03 (900023903| 
4] 0.4 | 0.09090205 | 0.49090205 | 1.49090205 0.147446 76 
5] 0.5 | 0.14744676 | 0.64744676 | 1.647446 76 0.220428 66 
6| 0.6 | 0.22042866 | 0.82042866 | 1.820428 66 0.311573 66 
7} 0.7 | 0.31157366 | 1.01157366 |2.01157366 0.422 788 89 
8] 0.8 | 0.42278889 | 1.22278889 | 2.22278889 0.55618171 
91 0.9 | 0.55618171 | 1.45618171 | 2.45618171 0.71408079 
10} 1.0 | 0.71408079 


26 


MST204 18.2 Tape 


Taulor Series method of order 3 | 


Compute the approximate solution bo the differential equation 


eae C 


usitin (0) =O 
method of order 3. 


0 
1 
2 
3 
4 
5 
6 
7 
8 
9 


10 


0.1 
0.2 
0.3 
0.4 
0.5 


0.6 


0.7 
0.8 
0.9 
1.0 


0.005 16667 
0.021393 35 
0.049843 22 
0.091 80173 
0.148689 54 
0.222076 72 
0.313 69845 
0.4254724 

0.55951791 
0.71817721 


eat <=O'l, 


ee 


yr = 


0.105 16667 


0.491 80173 
0.648 689 54 
0.822076 72 
1.013 698 45 
1.2254724 

1.45951791 


Oa. 


EN 


(9) 


ct Pe 
ie 


1.105 16667 


1.491 80173 
1.648689 54 
1.822076 72 
2.013 69845 
2.2254724 

2.45951791 


, |, using the Taylor Series 


0.005 16667 
0.021393 35 


0.148689 54 
0.22207672 
0.313698 45 
0.425472 4 

0.55951791 
0.71817721 
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We summarize the procedure for the Taylor series method of order 2 in the 
following procedure box. 


Taylor series method of order 2 
1. You are given a differential equation 


y’ = m(x, y) 


and a value for y(Xo). 


; . dm 
. Differentiate the expression for m(x, y) to obtain 7a y), 


remembering that y is a function of x. 


. To apply the Taylor series method of order 2, choose a step 
size h, and calculate Y;, Y>,... from 


Yai =¥+hy,’ + gh’y,” 
where 


Y,’ = m(x,, Y,) 
dm 

yy = -—(,), 

aoe: (Xr %) 


Yo = y(Xo) and Xp =Xo + rh. 


4. Y, is an approximation to y,. 


The procedure for using the Taylor series method of order n can be stated, in a 
similar way, as follows. 


Taylor series method of order n 
1. You are given a differential equation 


y’ = m(x, y) 
and a value for y(Xo). 
Differentiate the expression for m(x,y) n — 1 times to obtain 
ressions for a ) f diet ) = (x,y) 
ex 1 GT NY), Gor OV )o ++ 9 Grr OY) 
P a ae pent 
remembering that y is a function of x. 


To apply the Taylor series method of order n, choose a step 
size h, and calculate Y;, Y>, ... from 


2 


h h" 
/ " ba 3 ee 
Ki hthy ts w+ +K 


sa = m(X;, ¥,) 


y Pl in ee 
dx 


n-1 
yo — d"~*m 
. Pexaratied dx"7} 


¥5 = Vireo) and X, = Xo + rh. 


(x, ¥,) 


4. Y, 1s an approximation to y,. 


There is a decision to be made here whether to use a higher-order Taylor series 
method with a larger step size or a lower-order Taylor series method with smaller 
step size in order to achieve a required accuracy. Reducing the step size can reduce 
the error significantly, but requires correspondingly more computational effort to 


. dm 
The notation = (x,, ¥,) means 
= : 


d 
that we evaluate a (x, y) at 
x 


x =x, and y = Y,. 
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obtain approximations over a range of x-values. Adding an extra term in the 
Taylor series would probably not involve significantly more work provided that 
the derivative were easy to obtain. We shall say more about such comparisons in 
the next unit. 


One of the main drawbacks of the Taylor series methods is that they involve the 
use of higher derivatives which may be very complicated expressions. For example, 
you might like to convince yourself that the implementation of the Taylor series 
method of order 4 for the differential equation © 


y=/xr+y 
is by no means easy. (Do not spend more than a couple of minutes doing this.) 


In addition the work involved in inputting these higher derivatives into a 
computer or programmable calculator may be prohibitive. Consequently in the 
next unit we will look at other methods, which do not require these higher 
derivatives. 


Exercise 3 


Write down the recurrence relation for the Taylor series method of order 1. Does this 
method have another name? 


Exercise 4 _— 

Use the Taylor series method of order 2 to solve the differential equation 
y’ = 3y + sinx with y(0) = 0. 

Use h = 0.2 to compute approximate solutions at x = 0.2 and x = 0.4. 


Exercise 5 


Write down, the recurrence relation and expressions for Y,’ (in terms of x, and Y,), Y,” (in 
terms of x,, Y, and Y,’) and Y,” (in terms of x,, Y,, Y,) and Y,”) required to use the Taylor 
series method of order 3 to solve the differential equation 


y' =siny with y(0) = 1. 


Taking h = 0.1, compute the approximate solution at x = 0.1. 
[Solutions to Exercises 3—5 on p. 51] 


Summary of Section 2 


1. A root of the equation 


f (x) = @ 
is any value of x for which the equation is satisfied. 


2. The Newton-Raphson method, given an initial approximation Xp, close to the 
root of 


f(x) = 90, 
can be used to generate a sequence Xo, x;, X2, ... of successively better 
approximations to the root using the recurrence relation 

2 

r+1— ae ‘ 

TO) 
The error in x,,, is proportional to the square of the error in x, and so if the 
sequence converges it does so very rapidly. 


3. The Taylor series method of order 2 for solving first-order differential equations 
is given on page 27. 


4. The Taylor series method of order n for solving first-order differential equations 
is given on page 27. 


5. The selection of the appropriate order of the Taylor series method and the 
choice of step size are difficult to determine for a given differential equation. The 
drawback of the method is that higher derivatives may be complicated expressions 
which are tedious to work out and time-consuming to input to a computer or 
programmable calculator. 
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3. The catenary (Television Section) 


To illustrate how some of the techniques of this unit are used in practice, we are rope 

going to look at a particular problem which arises out of the study of catenaries. Ga ee 
A catenary is the shape taken by a freely hanging rope or cable fixed at both ends, 

as shown in Figure 1. Figure 1. A catenary. 


Here is a statement of the problem we are going to tackle: Take a length L of 
rope (where L is greater than 2m) and Suspend it between two points 2m apart in 
the same horizontal plane. What is the displacement d of the middle of the rope 
below this plane? (L and d are shown in F igure 2.) 


Figure 2. Television problem is to find d given L. 


This problem has applications in the design of telephone cables or electricity 
cables, in that temperature changes cause the cables to increase or decrease in 
length. 


It would be useful if we could find a formula expressing the displacement d in 
terms of the length L, because it would then be possible to determine d by 
substituting the given value of L in that formula. Unfortunately such a formula is 
impossible to find. However, we shall see that the Newton-Raphson method can 
be used to determine d for any given value of L. It is also possible to use Taylor 
polynomial approximations to determine an approximate relationship between d 
and L which gives good answers for small displacements. To tackle this problem 
we need two pieces of information: 


(i) the equation expressing the catenary as a function of x 
(ii) a formula for the length of a curve. 


These two topics are the subject of the next subsection, which you should read 
before watching the television programme. 


3.1 Hyperbolic functions 


The modelling of the catenary, which assumes that the mass per unit length is 
constant along the length of the rope, leads to a second-order differential equation. 
Unfortunately we do not have time to go into the details here. For a suitable They are givén, for example, 


choice of the x- and y-axes the equation of the catenary can be written as in Differential Equations with 
Applications and Historical 

1 - + =| Notes by G. F. Simmons, 

ee 


(1) McGraw-Hill, 1972. 
where a is a constant whose value is chosen to select the particular catenary we 
want from a whole family of catenaries. (For each length L of rope we will have a 
different value of a.) The suitable choice of axes (see Figure 3) is such that 


y'(0) =0 


, Qa 


1 
and y(0) = = 
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y(0) = 0 


=] 0 jee 
Figure 3. Axes for the catenary in the case a = 1. 


The expression in brackets in Equation (1) appears so often that it is given a 
name. We define the hyperbolic cosine of the number ax as 


cosh ax = eae 
2 
The corresponding hyperbolic sine is defined as 
em —_e &® 
a 
sinh ax 5 


The equation of the catenary can now be written as 
y= Pcie ax. (2) 


The hyperbolic functions have some very interesting properties, which we will use 
in the television programme. These properties are very similar to the properties of 
ordinary sines and cosines. In fact, to distinguish between the two types of 
function we call the well-known sines and cosines circular functions. This name 
arises because any point on a unit circle can be defined in terms of the cosine and 
sine of the angle 0 shown in Figure 4. 


Since the equation of the unit circle is given by 
xP ey? ee 
the well-known equation 
cos 6 + sin’ 6 = 1 
shows that the point (cos 0, sin 8) lies on this circle. 


To see why we call cosh and sinh hyperbolic functions, consider the curve, called a 
hyperbola, whose equation is 


x* — y* = |. (3) 


I shall show that the point whose coordinates are (cosh 6, sinh 6) lies on the 
hyperbola, i.e. that it satisfies Equation (3). The number @ can take any value, but 
its geometical significance is less obvious than for the circle. 


Now 
6 -6 6 4-6 
cosh 6 + sinh 8 = sc sett a oe = 
2 2 
and 
fae gn id 
bp sah} = bo ee ee 
cos sinh = (242°) - FS] we 
Thus 


cosh20 — sinh?9 = (cosh @ + sinh 6) (cosh @ — sinh 8) 
—  e — 


as required. 


The symbol cosh is 
pronounced as ‘cosh’. 


The symbol sinh is 
pronounced as ‘shine’. 


Figure 4. Any point on the 
unit circle can be expressed in 
terms of 0. 


(cosh@, sinh) 


Figure 5. Any point on the 
hyperbola x? — y? = 1 can be 
expressed in terms of 0. 
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Note the similarity between the result 
cosh7@ — sinh?6 = 1 

for hyperbolic functions and the result 
cos76 + sin?0 = 1 

for circular functions. 


The following table lists some properties of hyperbolic functions together with the 
corresponding properties of circular functions. They can be derived using the 
definitions of cosh x and sinh x in terms of exponential functions. If you have time 
you might verify some of these properties. You are not expected to memorize 
these; they will be included in the Handbook. 


Hyperbolic functions Circular functions 
e+e* e+e°* These first two expressions for 
cosh x = oe ae nieakecs bos asta “peices the circular functions were 
given in Unit 5, Subsection 
aoe ee ex 4.2. 
sinh x = ————_ sin xX = ——__—_— 
2 2i 


cosh?x — sinh?x = 1 cos*x + sin?x = 1 


2 x) = sinh x mane x) = —sin x 
dx dx 


x (sinh x) = cosh x fh (sin x) = cos x 
dx dx 


2 + 6 2 cs 6 8 


ee a = = er en ay = 

cosn xX = Ty gy te, Tt cos x = =e a a a 
. eee lee aoe <2 yo ee . 
sin 2 ae a Pe eS a ae 


The last two properties in the table can also be derived using Taylor’s Theorem, 
given in Subsection 1.3. For example, if f(x) = cosh x, the conditions of Taylor’s 
Theorem are satisfied for any interval [a, b] and we may write the 2nth Taylor 
expansion about 0, with remainder, as 


2 
cosh x = cosh 0 + x sinh 0 + 5 cosh 0 + -- 


x2" 2n+1 


a on 


where c,,. lies between 0 and x. 


sinh c, (4) 


Now cosh 0 = 1 and sinh 0 = 0. Hence Equation (4) simplifies to 


Se Cale = + oe 
SS a (2n)!  (2n+ 1)! 
As n becomes large we obtain the required expression for cosh x, since for any x 
the remainder term approaches zero for large n. (We will not prove this, as it 
would take us too far from the main subject of the unit.) 


sinh c,. 


In Exercise 3 you will be asked to obtain a formula for the length of a catenary. 
To do this you will need to use the following result. 


Length of a curve 
The length L of the curve y = f(x) between xo and x, is given by 


Xi 4 
L= | /1 - fa dx. 
os dx 


You can verify this result in 
the case of the straight line 

y = ax 
between xo = 0 and x, = 1 by 
showing that 


L=./1+a’. 
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For most curves a numerical method would have to be used to evaluate L, but in 
the case of a catenary we can determine L exactly. | 


Exercise 1 

Show that: 

(i) cosh (— x) = cosh x 
sinh (—x) = —sinh x 


- d 

(ii) Fy (oosh x) = sinh x 
, ae 
— (sinh x) = cosh x. 
dx | 


Compare these results with the corresponding properties of sines and cosines. 


Exercise 2 
Use Taylor’s Theorem to derive the Taylor series for sinh x as 
a 7 


4 h — = = ar 
sin er at a 
Exercise 3 
Determine the length of the catenary, between x = —1 and x = 1, given by 


1 
y =-—cosh ax. 
oo ae 


[Solutions to Exercises 1-3 on pp. 51-2] 


Now watch the television programme: ‘Catenaries — numerical approximation’. 


3.2 The rope problem 


TV 18 


As I stated in the introduction to Section 3, the problem considered in the 
television programme is to find the displacement d at the middle of a length L of 
rope suspended between points two metres apart as in Figure 2. 


In Subsection 3.1 I stated that the equation of a catenary is given by Equation (2) 
as 


1 
y= Pas ax. (5) 


-1 0 1 x 
Figure 6. The catenary. 


To obtain an expréssion for d we note that the height of the catenary above the x- 
; 1 
axis atx =lisd+ : (see Figure 6), ie. 
1 
yl) =d+-. 
a 
Using the equation for the catenary (5) we have 


y(1) = S onl a=d+ ss 
a a 
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Hence 


a= “(cosh a — 1). (6) 


This gives d in terms of a, whereas we really require an expression for d in terms 
of L. | 


From Exercise 3 we know that the length of the catenary is given by 
Bi 
Dos 7 sinh a. (7) 


Equations (6) and (7) define d and L in terms of a, but unfortunately there is no 
way to establish a direct formula for d in terms of L. It is possible to use these 
formulae, with the help of the Newton-Raphson method, although it does involve 
a lot of work. 


Example 1 
Determine the displacement d if L = 2.2 m. 


Solution 
From Equation (7) we have 


2.2 = * sinh a. (8) 


Rearranging this equation gives a as a root of the non-linear equation 
f (a) = sinha —1.1a=0. (9) 


From the sketch graph of f (a) in Figure 7 we can see that there are two non- 
negative roots for Equation (9). However, the root a = 0 is not a solution of 
Equation (8) and the value of a we require is approximately a = 0.76. 


Figure 7. Sketch of f(a) = sinh a — 1.la. 


The Newton-Raphson method, derived in Subsection 2.1, determines a root of the 
equation f(x) = 0 using the recurrence relation 


= iG) 
f'(,) 


Xp+1 = Xr 


with Xo given. 
For our function of a we have 
f(a) = sinh a — 1.1 a and  f’(a) = cosh a — 1.1. 


Hence the recurrence relation is 


oo sinh a, — 1.1 a, 
ee: cosh att | 
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With dao = 0.76 we have 
a, = 0.763 424 64 
a, = 0.763 400 80 
= a3. 


After only 3 iterations the sequence has converged to the required root. With 
a = 0.763 400 80 we can now use Equation (6) to evaluate d as 


1 
d = — (cosh a — 1) = 0.400601 6, 


Thus when L = 2.2m, the displacement at the centre of the rope is given as 
d = 0.401 m correct to 3 decimal places. 


In the above example we have managed to find d for a given length L. However, 
for each value of L we have to determine an initial value ao, use the Newton- 
Raphson method, and then substitute the computed value of a into Equation (6) 
to determine the displacement. 


In the second part of the television programme we look for a simple approximate 
relationship between d and L using Taylor polynomials, which works well for 
small displacements d. The advantage of a simple relationship is that it will be 
easier to make qualitative predictions about the behaviour of the rope for small 
displacements. Thus, although it is possible to obtain a more accurate 
approximate relationship between d and L, the simple one may be more useful in 
practice. 


From the properties of hyperbolic functions we know that 


2 + q® 


h mass 1 . ed oe 
cosh a = Ne 
Hence the 2nd Taylor polynomial about 0 for cosh a is 
cosh a ~ 1 + 4a’. Because the a® term is zero, 
| pe Paes this is actually the 3rd Taylor 
While a is reasonably close to zero this will provide a good approximation to polynomial as well as the 2nd. 


cosh a. In Equation (6) we have 


1 
d = —(cosh a — 1) 
a 


1 
~—-(1 + 4a? — 1) (using the 2nd (or 3rd) Taylor 
5 polynomial about 0 for cosh a) 
= $a. (10) 


Similarly we have 
ae 
sinha=ata tot 
and the 3rd Taylor polynomial about 0 for sinh a is 
sinh a ~ a+ ta’®. 


Thus the equation for L is 


Z 
L =-sinh a 
a 


2 
~—(a + ga’) 
a 
=2+4+ 4a’. (11) 


Now Equation (11) can be rearranged to give a in terms of L as 


a~ JL — 2). 


Hence using Equation (10) the approximate displacement d* is given by 


d* =4,/3(L — 2) = ./AGL — 6). (12) 
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Now we have an explicit expression for d* for any value of L. For example, if 
L = 2.2 we have 


d* = ./0.15 ~ 0.387, 
which is quite close to the more accurate value d ~ 0.401 m, obtained in Example 1. 


Comparing the approximate results, for different values of I, with the true results 
in Figure 8, we can see that for values of L slightly greater than 2 the 
approximation is very good, but as L increases the approximation gets worse. 


d 
; l 
true values Each of the true values was 
e computed using Equations (6) 
and (7) and the Newton- 
Raphson method. 
0.5 approximate solution 
0 


sh See 2.1 22 23 2.4 2:5 2.6 at i 


Figure 8. Graph comparing the approximate displacements with the true displacements. 


We conclude that our approximate method is a very useful way of determining the 
approximate displacement d* when d is small. Furthermore, the approximate 
relationship tells us something we might have otherwise missed. For L slightly 
greater than 2 there is a very sharp increase in the displacement. For example, if 

_ L= 2.001 we find that the displacement d ~ 0.027 m, i.e. if L is increased by 1mm 
from the zero-displacement position then the displacement will be 2.7cm, 27 times 
as big. Hence the problem of finding d given L is absolutely ill-conditioned for L 
close to 2m. The ill-conditioning is demonstrated at the end of the television 
programme, illustrating that approximate relationships can give us information 
that would otherwise be difficult to deduce. 


Exercise 4 


Use the Newton-Raphson method to determine the displacement d, to an accuracy of 6 
significant figures, when L = 2.08 m. 


Exercise 5 


Express a in terms of L if the 5th Taylor polynomial about a = 0 is used to approximate 
sinh a in Equation (7). 


Exercise 6 


Use the approximation derived in Exercise 5 together with Equation (6) to find the 
approximate displacement if L = 2.08 m. 


[Solutions to Exercises 4-6 on pp. 52-3] 


Summary of Section 3 


1. The equation modelling a freely hanging rope suspended between two points is 
2 | am + a 
oa 2 


1 
= -—cosh ax, 
a 


where a is a parameter which depends on (among other things) the length of the 
rope, and cosh is the hyperbolic cosine function. The curve with this equation is 
called a catenary. 
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2. Hyperbolic functions, given by 


2 
and 
e—e ~ 
Sse itis 
sinh x 5 


have many properties which are similar to those of the circular functions sin x and 
cos x. A list of these properties is given on page 31. 


3. The length of a curve between x» and x, is given by 


X14 Zz 
Lom 1+ fa dx. 
oS \ ae 
In most cases a numerical method would be required to evaluate this integral. 
4. The displacement d for a freely hanging rope of length L suspended between 
two points 2m apart is given by | 
i= Sod a—1) 
a 
where 


L= Saad a. 
a 


For a given value of L, a can be determined using the Newton-Raphson method 
and this value can be used to determine d. An approximate relationship between d 
and L, obtained using the 3rd Taylor polynomials about a = 0 for cosha and 

sinh a, is given by 


d* = AGL — 6). 


For 2 < L < 2.4 this gives approximations within 7 % of the true solution. 


4 Interpolation and integration 


4.1 Interpolation 


In Section 1 we derived the Taylor polynomials of degree n by assuming that we 
knew or could determine n derivatives of a function at a particular point «. 
However, in many problems the values of the derivatives are not known, but the 
data is given at a number of points Xo, X1, X2, X3, .-. , X». Diagrammatically we 
could represent this as in Figure 1. 


_ 9 x] Xx? X3 aa Xn x 
Figure 1. Data given at n + 1 points Xo, X1,... , Xr 


It is possible to show that there is just one polynomial of degree n which takes the 
desired value of f at each of the n + 1 points Xo, x1, ..., Xp. 
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Once we have determined this polynomial we can use it to approximate the 
function at other values of x. This process of determining the approximate value of 
the function at a point x between xo and x, is called interpolation. (The process of 
using the polynomial to determine approximate values outside this range is called 
extrapolation. Extrapolation is much less accurate than interpolation and we will 
not consider this topic here.) One method of computing the coefficients of the 
interpolating polynomial of degree n, p(x), given as 


P(X) = dg +ayx +... + a,x", 


is to set up a system of equations. Since p(x;) = f (x;) for i= 0, 1, ..., n, we can 
write 


P(Xo0) = Ao + AyXq +... + AnXq" =f (Xo) 
P(X1) = Ag + a4X, +... + a,x," =f (x) 


P(Xn) = Oo +-0,xX, +... + aX, = f (x,) 
This is just a system of linear equations, like those you met in Unit 9, Simultaneous 


linear algebraic equations, for the n + 1 unknowns do, a, ..., a,. In matrix form 
the problem is to find do, a;, ..., a, which satisfy the equations 


1 2 . Note that in this problem we 
a ay ao J (%o) know the values of xo, x, ..., 
1 X41 +" re be ee ay, tT (x1) Xn and f (Xo), f(x1),.-..f (%n) 
eae ; ; ee eas ee and we wish to determine 
° ; ° : Ap, A1,..-, Ay. 
Lo es le) 
Example 1 


Determine the interpolating quadratic polynomial which passes through the three 
points (0, 1), (1, 2) and (2, 4). 


Solution 
Let the interpolating polynomial be given by 
P(X) = dg + a,x + apx?. 
The three equations corresponding to the three data points are 
p(0) = ao + a,(0) + a2(0*) = 1 
p(1) = ap + a, (1) + a,(17) = 2 
p(2) = ap + a;(2) + a(2”) = 4 


giving 
Ao cae. 
Ag t a, + A= Z 
ag + 2a; ei 4a, = 4 
Le. 
t 0.9 ao l 
ee mee a,|= 2 
Pa die. (4 


Using Gaussian elimination, as described in Unit 9, we reduce this matrix to its 
upper triangular form as 


L- 0 Ola, 1 

i ee ele we ee 

8 2a. 1 
whence 


a, = 4, a, =4 and a = 1. 
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Thus the interpolating polynomial is 
p(x) = 1 + 4x + 4x?. 


Given a large set of points it is not always desirable to compute the interpolating 
polynomial of corresponding degree. This process can be very time-consuming, 
prone to the build-up of rounding errors in computation, and the resulting 
polynomial tedious to evaluate using a calculator. More frequently, given a value 
of x at which we need to approximate the function, we construct a polynomial 
using just a few of the data points. 


Example 2 
Here is a set of data obtained from four-figure tables for log, x. 


loge Xn 


Use an interpolating polynomial of degree 2 to obtain an approximation to 
log, (1.134). 


Solution 


For a quadratic polynomial we need to specify 3 points xo, x; and x2. We choose 
the three points closest to x = 1.134, namely 


Xo wy 4AZ, ¥, = FA. x, = 1.14. 
If p(x) = Ay + ax + a,x”, we obtain the following equations for do, a; and a2: 
p(1.12) = do + 1.12a, + (1.12)? a. = 0.1133 
p(1.13) = do + 1.13a, + (1.13)? a. = 0.1222 
p(1.14) = ao + 1.14a, + (1.14)? a. = 0.1310 


1 142 1.2544][ao] [0.1133 
1 1.13 1.2769 || a, | =| 0.1222). 
1 1.14 1.2996 |ja,} | 0.1310 
Using the computer package SIMLIN,I obtained the solution 
: ay = — 1.51635, a, = 2.015 08, a, = —0.50004. 
Hence 
p(x) = —1.516 35 + 2.01508 x — 0.50004 x’, 
giving 
p(1.134) = 0.1257 to 4 significant figures. 


It is worth pointing out that the problem of evaluating the coefficients of p in the 
above example is highly ill-conditioned. To illustrate this I repeated the 
calculations with the right-hand side values correct to five significant figures. The 
new set of equations to solve is 


1 112 1.2544][ao] [0.11333 
1 1.13 1.2769 || a, | =| 0.12222 |. 
1 114 1.2996 ||a,| | 0.13103 
SIMLIN this time gives 
dy = —1.38863, a, =1.78908 and a, = —0.40004, 


MST204 18.4 . 39 


so that 
p(x) = —1.388 63 + 1.78908 x — 0.40004 x?. 


This is completely different from the polynomial obtained in Example 2, and the 
difference is caused by very small changes in the data. The reason for this ill- 
conditioning is partly due to the fact that the rows of the matrix are almost 
linearly dependent (you can verify this if you wish!). However, if we use the new 
polynomial to evaluate an approximation to log, (1.134) we get 


p(1.134) = 0.12575 to 5 significant figures, 


and this answer is almost identical to the answer obtained in Example 2. The true 
solution (to 8 decimal places) is 


log. (1.134) = 0.125 751 21, 


indicating that the answers we obtained are remarkably good. Although the 
polynomials are completely different, the difference in their values for 

1.12 < x < 1.14 is less than 0.00004, ie. approximately the same as the differences 
in the data! 


What has happened is that although the problem of computing the coefficients of 
the interpolating polynomial is ill-conditioned, the problem of determining the 
approximate value for log, (1.134) is well-conditioned. This phenomenon occurs 
very frequently in interpolation problems where the x-values are close together. 
You need to be aware that this can occur, although it almost always makes no 
difference to the final answer. 


Exercise 1 


Determine the quadratic whose graph passes through the three points (— 1,2), (0,3) and 
(1,0). 


Exercise 2 
Given the data in Example 2, use an interpolating polynomial of degree 1 to compute an 
approximation to log, (1.134). 


[Solutions to Exercises 1 and 2 on p. 53] 


In Section 1, when we discussed Taylor polynomials, we stated Taylor’s Theorem 
which could be used to determine the error in the approximation. In the same way 
there is a very useful theorem which can be used to determine the error in the 
interpolating polynomial. In this theorem we assume that all arithmetic has been 
carried out exactly and there are no rounding errors in the values 


I Bink 7 akc, ). 


Theorem 4.1 
You are given a function f(x) and an interpolating polynomial 
p(x) of degree < n such that 


P(x;) = f (x;) for i= 0, 1, 22+ A, 


where the numbers xo, x, ..., X, are arranged in increasing order. 
Then the error, e(x) = p(x) — f(x), where x lies between x, and x,, 
is given by 


ie 
(n+ 1)! 


where c, is some number satisfying x9 < c, < x, 


(1) 


e(x) = —(x — Xo)(x — x1) +++ (x — x) 


You can see that the presence of the product (x — x9) (x — x;)...(x — X,) ensures that 


&(x;) = 0 int i= 0, 1. 8. 
as required, while the term 
RES 
(n+ 1)! 


is very similar to the corresponding term in Taylor’s Theorem. 
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The following procedure uses Theorem 1 to compute an error bound at a particular 
value of x between xX, and x,,. 


Error bounds for interpolating polynomials 


1. Let p(x) be the polynomial of degree <n which interpolates the 
function f(x) at Xo, X1,..., X,, Where x9 < X1 <...< X,,. 


2. The error é(x) = p(x) — f(x) is given by 
theese 2 


&(x) a — (x a Xo) (x a Ma) AX aig 2) DE 


where c, lies between Xo and x, 
3. Determine a value for M such that 
SPO <M 
for all values of c between x, and x,,. | 
4. A bound on the error e(x) is then 


\Gx — x0)(x — x1)...(% — x0) M 


le(x)| < ar 


To see how this procedure can be used, we calculate an error bound in the 
following example. 


Example 3 
(i) Show that |1/c?| < 1 for all values of c between 1 and 2. 
(ii) The linear approximation to log, x whose graph passes through the points 


(1, 0) and (2, 0.6931) is given by 
p(x) = 0.6931x — 0.6931. 


Determine an error bound at x = 1.5 and compare this bound with the 
actual error at x = 1.5. 


Solution 

(i) Since c > 1, we have 0 < 1/c < 1, and so |1/c?| < 1. 

(ii) Using the procedure for computing error bounds, we note that with 
Xo = 1 and x, = 2 the error function is given by 


e(x) = — Ge — 1) — 29 


where c,, lies between 1 and 2. 
Now f(x) = log, x, and so 
f'(x)=1/x and f"(x) = —1/x’. 


Hence from part (i) an upper bound for | f”(c)| = |1/c?| where 
1 <c <2 is given by 


M = I. 
Thus |e(x)| < aes sao and at x = 1.5 we have 
le(1.5)| < IG.5 = DG5 — 2) == 195. 


Z 


Thus an error bound for this approximation is given by 0.125. The actual 3 
error is given by 


p(i.5) — log, 1.5 = 0.34655 — 0.405 47 
= —0.058 92 


so that |e(1.5)| = 0.058 92. This is certainly less than the error bound 0.125. 


The only error in the 
calculation of p(x) is in giving 
log, 2 to only 4 decimal 
places. - 
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Exercise 3 

(i) Show that, for —1 <c < 1, we have |e| < e = 2.7183.... 

(11) The function f(x) = e* is to be approximated by the interpolating quadratic 
polynomial which interpolates f at x» = —1, x; = 0 and x, = 1. Determine a 


bound for the error in this quadratic at x = 0.5, ignoring rounding errors. (Note 
that you are not asked to compute the interpolating quadratic polynomial.) 


[Solution on p. 53] 


4.2 Methods of numerical integration 


In the last subsection we looked at methods of interpolation. In this final 
subsection we examine how these interpolation methods can be used to derive 
approximations to the integral 


f= [ food 


where f is some given function. 


By integrating the interpolating polynomial, which is a fairly straightforward 
process, we can obtain formulae for the approximation of the integral. Such 
formulae are known as Newton-Cotes formulae. The methods derived here, Euler’s 
integration method, the trapezoidal method and Simpson’s method, will be used in 
the next unit to solve differential equations, but they are clearly useful in their own 
right as well. 


Euler’s method 
Suppose we want to approximate the integral 


I : I (x)dx. 


Figure 2 shows one simple approximation: we replace the function f (x) by the 
constant f (xo) for all values of x between xo and x,. This is equivalent to using a 
very crude ‘approximating polynomial’, given by 


P(x) =f (Xo). 
The integral 


r= [ feoae 


is the area under the curve y = f(x) between x, and x,. We approximate this area 
by the area under the approximating polynomial as 


1~ { | peeks { 1 i Goldx = @: conor 


x 


= Af (Xo) where h = x; — Xo, 


i.e. the area under the curve y = f(x) is approximated by the area of the rectangle 
of height f(x ) and width h. | 


Clearly this is a very crude approximation, and Euler’s method would rarely be 
used in this form. However, we can get more accurate results by using a composite 
form of the method. Suppose we wish to approximate the integral of f over the 
interval [a, b]. Then we divide this interval into n subintervals of equal width h 
and use Euler’s method on each subinterval (see Figure 3). 


The approximation to the integral is then 


x 


b x4 n 
[ foods = | SI (x)dx + +++ + tT (x)dx 


~ hf (Xo) + hf (x1) + +++ + hf (xq-1) 
=h(f (x0) +f(x1) + °° + f (%q—1)). (2) 


The reason for calling it 
Euler’s method is that it is 
Closely related to Euler’s 
method for differential 
equations. See Unit 19. 


ba 


function y = f(x) 
approxi 
polyno 
y = pt 


Xo xX] x 


Figure 2. Euler’s method. 
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a=Xv x) x9 abe rh eee xX, =5 »¢ 


Figure 3. Euler’s method in composite form. 


Example 4 
Determine the approximate value of 


2 
1= | log, x dx 
1 


~ using the composite Euler’s method with h = 0.1. 
Solution 
Using Equation (2) with h = 0.1, we have 


I ~0.1 (dlog.1 + log.1.1 + log,1.2 + --- + log.1.9) 
= (0.351 22. 


For this example the true solution is given by 


2 
log. x dx = [x log.x — x]? 
1 
= 2log,2-—2+1 
= 0.386 29. 
The approximation is not particularly good, but from Figure 3 we can deduce that 


if we use a very small interval width h we can achieve much better accuracy. 


Euler’s method in composite form 
1. You are asked to obtain a value of the integral 


f= f foodx 


2. Choose values for h and n such that 
b=a-+nh. 

3. Approximate J using Euler’s method in composite form as 
I = Wf) +f) 42-9 OE 


where x, =a-+rh. 


Exercise 4 


Use Euler’s method in composite form with h = 0.2 to approximate the integral of e* 
between x = 0 and x = 1. 


[Solution on p. 53 | 


The trapezoidal method 


Consider the linear interpolating polynomial approximation to the function 
y =f (x) which passes through the data points (xo, f (xo)) and (x1, f(x1)). 


The integral 


ia [ Fix 
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is approximated by the area under the graph of the interpolating polynomial. For 
linear interpolation this area is a trapezium—a quadrilateral with two opposite 
sides parallel. The area of this trapezium is the product of the width h = x, — xo 
times the average height, which is 4(f (xo) + f (x,)). The approximation to the 
integral J is thus 


I = $h(f (x0) +f (x1)). (3) 


function y = f(x) 


interpolating 
polynomial 


X90 
Figure 4. The trapezoidal method. 


This method of approximating the integral is known as the trapezoidal method. We 
would expect the trapezoidal method to give more accurate results than Euler’s 
method, since the trapezium fits the area under the curve better than the rectangle 
does. 


Example 5 y 


; : y = log.x 
Determine an approximate value of 


z 
i= | log, x dx 
1 


using the trapezoidal method. 

Solution 

Equation (3) gives the approximation as Figure 5. The trapezoidal 
I ~ #(log, 1 + log, 2) = 0.34657. method tor Example'5. 


The correct solution is given in Example 4 as I = 0.386 29. Even with this fairly 
crude method we have obtained a reasonable approximation to the integral. 


The trapezoidal method is most commonly used in composite form. That is, if we 
wish to approximate the integrand f over the interval [a, b], then we divide this 
interval into n subintervals of equal width and use the trapezoidal method on each 
subinterval, as shown in Figure 6. 


function y = f(x) 


a=Xo xX] xX) mee 
Figure 6. The composite trapezoidal method. 


= | fxd ~ SoC (0) + f0%)) + HAC (a1) + fl) + 
+ xh io) + f (xn) 
= (SF (0) +f (%1) +f(%2) + +f (%n—1) + 4£ en) (4) 
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Example 6 
Determine the approximate value of 


02 
l= log, x dx 
1 


using the composite trapezoidal method with h = 0.25. 
Solution | 
Equation (4) with f(x) = log.x and h = 0.25 gives 


I ~ 0.25 (log, 1 + log, 1.25 + log, 1.5 + log, 1.75 + $log, 2) 
= 0.383 70. 


This is almost correct to 2 decimal places, since the true solution is J = 0.386 29. 


The trapezoidal method in composite form 


1. You are asked to find a value of 


l= { f (x)dx. 


2. Choose values for h and n such that 


b=a-+nh. 


3. Approximate J using the trapezoidal method in composite 
form as 


I ~ h(ES (x0) +f (x1) + f(%2) + +L %n-1) + EF %n)) (5) 


where x, = a + rh. 


Exercise 5 


Use the composite trapezoidal method (Equation (5)) with h = 0.2 to estimate the integral 
of e* between x = 0 and x = 1. Compare your answer with the true solution and the 
solution obtained in Exercise 4. 


[Solution on p. 53] 


Simpson’s method 


Our third method of numerical integration is known as Simpson’s method. It 1s 
derived by interpolating three equally spaced data points (xo, f(xo)), (x1, f(*1)) 
and (x2, f(x2)) by a quadratic polynomial and then integrating this quadratic, as 
in Figure 8. Since it can allow for the curvature of the graph of the function f, we 
could expect Simpson’s method to give much more accurate results than Euler’s 
method or the trapezoidal method. 


y 


function y = f(x) 


interpolating 
quadratic polynomial 


Figure 8. Simpson’s method. 


To simplify the derivation of the formula for Simpson’s method we make a change 
of variable, writing 


s=x—%X, | (6) 


N 
Figure 7. The composite 


trapezoidal method applied in 


Example 6. 
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and we denote the quadratic interpolating polynomial, written in terms of the 
variable s, by 


q(s) =a + bs + cs’. 


From Equation (6) and the conditions determining the interpolating polynomial 
we see that: 


when x = Xo we have s = —h and q(—h) =f (Xo), 
when x = x, we have s = 0 and q(0) =f (x,), 
when x = x2 we have s = h and q(h) = f (x2). 


Hence the conditions determining the coefficients a, b, c in the polynomial q(s) are 


q(—h) =a — bh + ch? =f (xo) : (7) 
q(0)=a = f (x1) (8) 
q(h) =a + bh + ch? =f (x,). (9) 


Equation (8) gives a = f(x,). Thus Equations (7) and (9) can be rewritten as 


—bh + ch* = f (xo) —f (1) 
bh + ch* =f (x2) —f (x1). 


These two equations can now be solved for b and c, giving 


» £2) —F (0) 
2h 


and 
£2) = 2 Fes) + Fo) 
= 2h? 
Thus 


f (x2) —f (Xo) 


f (x2) —2 f(x1) +f (Xo) . 
2h 


2h?. 


S+ 


q(s) =f (x1) + | 
To approximate the integral of f(x) between x, and x, we write 


) -Sesydx ~ | plxax 


= [aod 


ss f as fea rea, 


1 (fea - 24 (%1) +i) = = 
2h 
ss | Fees " [foe fel . 
, (£@2) = 2F 1) +F60)) 8°) 
2h? Stee 


= hf (1) +0 +5 (flea) — 2F le) + f%0)) 


Hence 
x2 h 
[ seedx = 5 (Fea) + 4f Ces) +f (o)). (10) 


Now that we have the formula for Simpson’s method, given by Equation (10), we 
can simply use this to approximate integrals. 


Do not spend a lot of time 
working through the 
derivation of Equation (10). 


When x = Xo, s = —h, and 
when x = x2, s = h. Equation 
(6) gives ds = dx. 


46 MST204 18.4 


Example 7 
Determine an approximation to the integral 


2 
i= | log, x dx 
1 


using Simpson’s method. 


Solution 
To use Simpson’s method we need two equal intervals of width h = 3. From We have not drawn a figure in 
Equation (10) we have, since h = 4(x2 — Xo) = 4, this case, as there would be no 
visible difference between the 
I ~ ¢ (log.1 + 4log.1.5 + log.2) graph of y = log, x and the 
— 1.385 83 interpolating quadratic. 


and this is correct to three decimal places. 


Like Euler’s method and the trapezoidal method, Simpson’s method can be used 
in composite form. We divide the interval [a,b] into an even number of 
subintervals and approximate the integral as the area under the quadratic for each 
pair of subintervals. 


y function 


Figure 9. The composite Simpson’s method. 


: h 
{ f (x)dx ~ 5 (F(0) + 4f(x;) + f(%2)) + 3 2) + 4f (x3) +f(x4)) + °°: 
h 
5 (f (x,-2) ae 4f (Xn-1) +f (Xn)) 


h 
= 3 Oo) + 4f (x1) + 2f(%2) + 4h (x3) + 2f(%4) + °° 
+ 4f (Xn—1) +f%n)). (11) 


The composite Simpson’s method is one of the most popular methods for 
numerical integration. We summarize the procedure in the following box. 


Simpson’s method in composite form 
1. You are asked to obtain a value for the integral 


P= { f (x)dx. 


2. Choose values for h and n such that n is even and 
b=a-+nh. 


3. Approximate J using Simpson’s method in composite form as 


1 = (F%0) + Af Cex) + 2F Cea) + 4f ts) + 2f xa) + 
Es 4f (Xn-1) +f (Xn) 


where x, = a+ rh. 
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Example 8 
Determine an approximate value for the integral 


y 
i= | log.x dx 
1 


using the composite Simpson’s method with h = 1. 
Solution 
With h = 4, Equation (11) gives 


I = yz (log-1 + 41og,1.25 + 2log.1.5 + 4log,1.75 + log,2) 
= 0.386 260 


and this solution is correct to four decimal places. 


Exercise 6 
(1) Use Simpson’s method to obtain an approximation to the integral of e* between x = 0 
and x = I. 


(i!) Use the composite Simpson’s method for the same integral with h = § and h = $. 
Check the accuracy of your solutions. 


Exercise 7 


The length of the curve y = sin x between 0 and : is given by 


m/2 
L= | (1 + cos? x)!/? dx. See Subsection 3.1. 
0 


Use the composite Simpson’s method with h = = to find the approximate length. 


[Solutions to Exercises 6 and 7 on p. 53] 


Summary of Section 4 


1. Given a set of data (xo, f(xo)), (1, f(X1)), --- 5 (Xm S(%n)), interpolation refers to 
any method of estimating the value of f(x) at an intermediate point x by 
constructing an interpolating function whose graph passes through these 
points and then evaluating this function at x. 

2. To construct an interpolating polynomial p of degree <n given by 

P(X) = dg + ayX + anx? +--+ + .,x" 


which passes through the points (xo, f(Xo)), (X1,f(%1)),..-5 (Xn, f(x")), we solve 
the n + 1 linear equations for ao,...,a, given by 


Mi) = J) ba 


3. The interpolation error for a polynomial of degree n is given in Theorem 4.1. A 
bound for this error can be found using the procedure on page 40. 


4. KEuler’s method for approximating the integral of f between x, and x, is given 
by 


| Fk = afin 


Euler’s method in composite form for approximating the integral of f between 
Xo and x, is given by 


| Nee iho) +f) Lf) + 4 Fe) 


where h = X;4; — X; i=... 2 — I. 
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5. The trapezoidal method for approximating the integral of f between xo and x, 
is given by 


[ feede = 5 (F060) $F) 


where h = x; — Xo. 


The composite trapezoidal method for approximating the integral of f between 
Xo and x, is given by 


™ fiosydx 2 OF Oo) + f0e1) + faye 47 D4 FED) 
where h = X;41 — Xi, i=0,....n—1. 


6. Simpson’s method for approximating the integral of f between xo and x, is 
given by | 


[Ferd ~3 (Fe00) + Af) +f) 


where h = x, — X; = X1 — Xo. 
The composite Simpson’s method for approximating the integral of f between xo 
and x, (n even) is given by | 
Xn h 
[ “seed ~ 3 (F(%0) + Af xs) + 2F (2) + Af rs) + + 2F Cea) 
+4 f(x,-1) +f) 
where h = X;11 — X;, | eee ee 2 


Simpson’s method is usually much more accurate than the trapezoidal method or 
Euler’s method. 


5 End of unit exercises 


Exercise 1 
(i) | Show that the maximum value of 2 tan csec’c is 0.20 to 2 decimal places for Hint: 
0<c<01. sec?x = 1 + tan?x 


(ii) | Determine the linear Taylor polynomial about 0 for tan x and give an error bound 


d 
oe —-tan x = sec*x 
on the approximation of tan x by this polynomial for 0 < x < 0.1. dx 


Exercise 2 — xX = sec x tan x. 
Determine the two positive roots of the equation 
f(x) =e& — 3x =0 


to 5 decimal places. 


Exercise 3 


Use the Taylor series method of order 2 with h = 0.1 to obtain approximations to the 
solution of the differential equation 


y=xy+1 with y(0) = 1 


at x = 0.1 and x = 0.2. Give your approximations to 3 decimal places. 
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Exercise 4 
(i) You are given just the following information about a function ge 
f (1) =1 and f (24) = 14. 
Use an interpolating polynomial of degree one to approximate the function at 
x = 1.44. 
(ii) | Verify that for 1 < c < 24 we have [4c~9/7| < 4. 
(iii) The function fin this problem is actually given by 


Save Jf % 


Determine an error bound for the approximation at x = 1.44 and compare this 
bound with the actual error. 


Exercise 5 
Use Simpson’s method in composite form with h = 0.25 to approximate the integral of 
tan x between 0 and 1. 


Exercise 6 ‘ 


Use Simpson’s méthod in composite form with h = 0.25 to determine the approximate 
length of the curve y = e* between x = 0 and x = 1. 


[Solutions to Exercises 1-6 on pp. 53-5] 
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Appendix: Solutions to the exercises 


Solutions to the exercises in Section 1 


1. (i) Differentiating term by term, we have 
p’(x) = 8x? + 9x? + 2x — 7. 
(ii) Integrating term by term gives 
[pcode = ax > sy + <— — — +2x+C 
where C is the constant of integration. 
(iii) Using nested multiplication, we have 
tig = 2 
u, = 0.29u9 + 3 = 3.58 
uz, = 0.29u, + 1 = 2.0382 
uz; = 0.29u, — 7 = —6.408 922 
ug = 0.29u3 + 2 = 0.141 412 62. 
Hence p(0.29) = 0.141 412 62. 
(iv) From (i), p’(x) = 8x? + 9x? + 2x — 7. 3 
Using nested multiplication for this cubic polynomial gives 
i= 3 
u, = 0.77uy + 9 = 15.16 
u, = 0.77u, + 2 = 13.6732 
u3 = 0.77uz — 7 = 3.528 364. 
Hence p’(0.77) = 3.528 364. 
2. Since f(x) = (1 + x)t2, we have 
fo) =H +x)? 
fx) = &(-)0 + x) 7? = -40 + x)”. 
When x = 0, (1 + x} = 1 for any k. Thus 
f (0) =1,f'(0) =4 and f"(0) = —4. 
~ Hence the 2nd Taylor polynomial is 
p(x) = 1+ 4x — $x’. 
At x = 0.44 we evaluate p(x) using nested multiplication as 
Up = —0.125 
u, = 0.44u,y + 0.5 = 0.445 
u, = 0.44u, + 1 = 1.1958, 
i.e. p(0.44) = 1.1958. 
This is correct to two decimal places, since 
(1 + 0.44)!/* = 1.2. 
3. f(x) =log,(1 + x). Differentiating gives 
: 1 1 
f(x) = oe f"(x) = me 
ed)! 


my ane 2 ie Ne fae n+ 
P= as POH oe 


(1 + x) 
At x = 0 we have 

f()=log.1=0, fO=1, f"O= —1, 

{ Ths 2 {0) =(-1"* (0 - 1)! 
Hence the nth Taylor polynomial about 0 for log,x is 
t=1)x" = 2x" 

2! 3! : n! 
ee nt t 5 
So 
n 


=x —4x?+47—--- 


as required. 


= (—1)"*! (n—1)!x" 


4. (i) Since 0 <c < 0.5, we have 1 <1 +c < 1.5, and hence 
ll +cl < 1.5. | 

(ii) f(x) =(1 +x). Thus 
f'(x)=S5+x)*, f(x) = 200 + x)’, 
f(x) = 6011 +x)? and f(x) = 120(1 + x). 


The error function for the 3rd Taylor polynomial is given by 
x* : 
e(x) = nay oe (cx) ag —5x4(1 + Cx) 


where c, lies between 0 and x. 


Now we are given that 0 < x < 0.5, so it follows that 
0 <c, < 0.5 and hence (by part (i)) that |1 + c,| < 1.5. Using 
this result to calculate a bound for «e(x), we obtain 


le(x) = 5x* SS 75x* Oc x= 05) 

which is the required result. 

5. (i) We know that —1 <sin c < 1 for all c, and hence 
lsin c| < 1. 

(ii) Let f(x) = sin x. Differentiating gives 


f(x) =cosx, f"(x)= —sinx, 


f(x) = —cosx, f(x) =sin x, 
f(x) =cosx and f(x) = —sinx. 
At x = 0, sin x = 0 and cosx = 1. Hence 
fO=9,, (O=1 FO=h 
f"O)=-1, f%O)=0 and f(O)=1. 
Hence the 5th Taylor polynomial about 0 for sin x is 
eee 
p(x) =x — 31 + 51 


The error is given using Taylor’s Theorem, as 


x 


s(x) = an Cc where c, lies between 0 and x. 


From part (i) we know that |sin c,| < 1. It follows that 
6 
x 
le(x)| < Sl 
Thus, at x = 1, we have 


1 
le(1)| < 6! = 0.001 388 89. 
Our error bound for the Sth Taylor polynomial about 0 
approximating sin x at x = 1 is thus 0.001 39, to three 
significant figures. (There is no point in quoting errors and 
error bounds to a large number of significant figures.) 


The actual value of sin 1 is 0.841 47099, and of the Taylor 
polynomial is 0.841 666 67, so the actual error is only 
0.000 196 (again to three figures). 


Solutions to the exercises in Section 2 


1. From the graph the three roots are approximately —0.5, 
—3.5 and —9.5. Using these values as initial estimates for 
the Newton-Raphson method with the recurrence relation 


(3x, + 27)x, + 40 


Xrt+i1 =X, — 
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we have the following table of results. (Note that the 
recurrence relation has been written in a form suitable for 
nested multiplication.) 


%» (Xo'= —0.5) | 269 = ~3.5) 


0; —0.5 ~~ 35 Oo 

1 | —0.497064 22 — 3.546 760 6 — 9.457 0507 
2 | —0.497068 01 — 3.546 402 3 — 9.456 529 8 
3 | —0.497068 01 — 3.546 402 3 — 9.456 5297 


As you can see, after only 3 iterations the results agree to 7 
or 8 significant figures. Hence, to 5 significant figures, the 
roots are —0.49707, — 3.5464 and —9.4565. 


2. (i) The equation 
xe *—0.1=0 


can be multiplied by 10e*, without introducing any new 
roots, to give 


10x —e* =0. 


0 
0.258 669 33 


0.02 
0.106 855 36 


Thus the original equation has roots whenever 10x = e*. 
Hence if we sketch the graphs of y = 10x and y = e* the 
intersections of the two graphs will give the roots (see Figure 
1). 


y) 


40 
30 


rook | 


Figure 1. Sketch of y = 10x and y= &. 
There are two roots at approximately 0.2 and 3.5. 


(ii) The Newton-Raphson method for this problem uses the 
recurrence relation 


x, exp (—x,) — 0.1 


7 feet an ay 


ie. 
x, — 0.1 exp (x,) 
Xp+1 =X, — ——_—_.. 
1 — x, 

RS ree | RANG = 39) 
0} 02 ao 
1 | 0.102675 34 3.575 3819 
2 | 0.111 74432 3507 1544 
3 | 041193255 3377 1352 1 
4 | 0.111 83256 3:57) 1521 


To 6 decimal places the roots are 0.111 833 and 3.577 152. 
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3. The recurrence relation for the Taylor series method of 
order | is given by 


Y+1= AB ua 
This is Euler’s method (Unit 2, Subsection 22). 


4. The recurrence relation for the Taylor series method of 
order 2 is given by 


Yei=¥+h¥ +4hh’y,” 
Now the differential equation is 
| y = 3y + sinx. 
Differentiating gives 

y” = 3y’ + cosx. 
Thus we have 


Y,’ = 3Y, + sinx,, 
Y,” = 3Y,’ + cosx,. 


Y,” = 3Y,’ + cos x, 


1 
1.756 0746 


0.02 
0.106 855 36 


The approximate solutions at x, = 0.2 and x, = 0.4 are 0.02 
and 0.106 855 36 respectively. 


5. The recurrence relation for the Taylor series method of 
order 3 is given by 


Yi =Y+hy + hn? y” +40 ¥” 
The differential equation is 
y' =sin y. 
Differentiating twice gives 
y" = y'cos y, 
y” =y" cosy — y’siny. 
Hence the required equations are 
Yo1=Y¥+hyY' + 5h’ Y” + th? y.”, where 


y= 1 
¥ =sin Y. 
se — Y,’ cos Y, 


y” = Y," cos ¥, = ¥.* sin ¥: 
For Xo = 0, Yo = 1 and h = 0.1 we have 
Yo = sin Yo = 0.841 47099 
Yo = Yocos Yo = 0.841 47099 x cos1 = 0.454648 72 
Yo” = Yo’ cos Yo — Yo" sin Yo 
(0.454 648 72 x cos 1) — (0.841 47099? x sin 1) 
= —0.350175 49. 


Thus 


Y, = 1+ (0.1 x 0.841 47099) + (0.005 x 0.454 648 72) 
— ( x 0.001 x 0.350175 49) 


= 1.086 362. 


Solutions to the exercises in Section 3 


e+e ™* 


1. (i) coshx = . Thus, 
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—x e 
cosh (—x) = = = = cosh x. 
—e* 
sinh x = ae Thus, 
ge: 
sinh (— x) = 5 = —sinh x. 


The corresponding results for sines and cosines are identical, 
i.e. 


cos (— x) = cosx and sin(—x) = —sin x. 
oe d je~+e™* e—e* 
(ii) Fx Osh x) = x | 5 | ee ee sinh x, 
’ d je*—e™* e+e* 7 
ax (sinh x) = re ae gee eS 4 cosh x. 


The corresponding results for sines and cosines are almost 
identical, i.e. 


= (sin x) = cos x and he (cos x) = —sin x. 
dx dx 


2. Taylor’s Theorem for f(x) = sinh x gives the 2nth Taylor 
ae about 0, with remainder, as 


x? 
sinh x = sinh + xcoshO + —sinhO += 31 Sake 4 - 
x2" 2n+1 
h —— cosh 
+ On mia 0+ (n+ a c. 


where c, lies between 0 and x. (Use the result of Exercise 1 
to do the differentiation.) With sinh 0 = 0 and cosh 0 = 1 we 
have 


3 x? x2n-1 
sinhx =x +37 carrey " Ga it 
x2nti 
h 
— 


As n becomes large we obtain the required expansion for 
sinh x, assuming (as is in fact true) that for fixed x the 
remainder term approaches zero for large n. 


1 
3 y= 7 cosh ax. Thus, 


d 
psf = sinh ax. 


dx 


Using the formula for the length of a curve, given in the box 
on p. 31: 


L= a A + ral dx 
= [ 1 + sinh? ax dx 
-1 
A 5 
= | ~ cosh ax dx (since cosh? ax — sinh? ax = 1) 
=i 
1 1 
= F sinh ax| 
a -1 
ie 
oe (sinh a — sinh (—a)) 
o. 
mae sinh a (using the result of Exercise 1). 


4. Substituting L = 2.08 into Equation (7) gives 


2.08 = * sinh a, 
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which can be rearranged as 
sinh a — 1.04a = 0. 
y 


y= Sinha -|-Oha 


Figure 2. Sketch of y = sinha — 1.04a. 


From the sketch graph of y = sinha — 1.04a in Figure 2, we 
deduce that there is a root at approximately 0.49. The Newton- 
Raphson method gives 
eeee sinh a, — 1.04a, 

Siete: cosha, — 1.04 ] 
ao = 0.49 
a, = 0.487 030 33 
az = 0.487 002 63 = a3. 


Hence we have a = 048700083. 


From Equation (6) we have 
1 
d = —(cosha — 1) 
a 


= 0.248 352 16. 
Thus the displacement for L = 2.08 m is, to 6 significant figures: 
d = 0.248 352 m. 
5. The 5th Taylor polynomial for sinh a about 0 is 
Rog 


p(x) =a+t t 790 


Now to obtain a in terms of L in this approximation, we note 
that 


2 a* 
37 0 
is a quadratic in a’, ie. 

oe 

+ ea 

ee aig 


The formula method gives 


a ee 
-3t /s+o-2) 
+ See 
: 1/30 


—10 + ./100 + 60L — 120 
= —10+ ./60L — 20. 


As a? must be positive there is only one possibility, i.e. 
Pp 


a* = —10+ 2,/15L —5 


giving 


a=./-10+2/ ISL — 5. 
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6. If ZL = 2.08 m we have from Exercise 5, remembering that it | 


is an approximation: 
a = 0.487018 59. 
Equation (6) gives 


1 
d = —(cosh a — 1) 
a 
~ 0.248 360 61. 


Comparing this with the accurate result in Exercise 4, we see that 
the approximation is correct to 4 significant figures, and even the 
fifth significant figure is nearly correct. 


Solutions to the exercises in Section 4 
1. Let the equation of the quadratic be 
p(x) = dy + 44x + ax?. 
Since p(x) must pass through the three data points, we must have 
p(—1) =a) — a, +a, =2 
p(0) =a 2 
P(l) =ajo+a,+a,=0. 


We can see that the solution is ag = 3,a,; = —landa, = —2,so 
that 


p(x) = 3 — x — 2x?. 
2. The interpolating polynomial is of the form 
P(x) = ap + a,x. 


The closest data points to x = 1.134 are x = 1.13 and x = 1.14. 
Thus we require p(x) to satisfy the equations 


p(1.13) = ap + 1.13a, = 0.1222 
p(1.14) = ap + 1.14a, = 0.1310. 


This gives a, = 0.88 and ay = —0.8722. Hence 
p(x) = —0.8722 + 0.88x 

and 
p(1.134) = 0.125 72. 


3. (i) The maximum value of |e‘| in the interval -1<c<1 
occurs at c = | since e is a positive, increasing function in this 


interval. Hence in this interval, 
0 < |e| < e = 2.7183 to 4 decimal places. 


(ii) The equation for the error is given by 


g(x) = —(x — xXo)(x — x1)(x — x5) = 


where c,. lies between xy and x3. 
Now f (x) = e*. Therefore f’’(x) = e*. Hence, 


a ink = 1) SPI 


Now for —1 <c, < 1 we have |exp (c,)| < 2.7183, and so 
2.7183 
le(x)| < |(x + 1)x(x — re 
For x = 0.5 this gives 


2.7183 
|e(0.5)| < |1.5 x 0.5 x TBs eer ee 


= 0.170 


Thus an error bound for the interpolating polynomial 
approximation at x = 0.5 is 0.170. 


to 3 decimal places. 
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4. Euler’s method in composite form for the integral of e 
between x = 0 and x = 1 with h = 0.2 gives 


1 
{ edx = 0.2 (e° + e%? + e%4 4 06 4 90-8) 
0 


= 1.5522 
(The true solution is e — 1 = 1.7183 to 4 decimal places.) 


to 4 decimal places. 


5. The composite trapezoidal method for the integral of e* 
between x = 0 and x = 1 with h = 0.2 gives 
1 
{ edx ~ 0.2 (5 e° +°e9-2 ake e?-4 a e?-6 4. e98 a te) 
0 


= 1.7240 


The true solution is e — 1 = 1.7183 to 4 decimal places. The 
result is much more accurate than the result obtained using the 
composite Euler’s method. 


to 4 decimal places. 


6. (i) To use Simpson’s method to approximate the integral of 
e* between x = 0 and x = 1 we have h = 0.5, giving 


0.5 6 p05... UA 
[ ear~Pee + 4e"° + e*) 
0 


= 1.7189 


This is much more accurate, for much less work, than the 
approximations in Exercises 4 and 5. 


(ii) With h = 4 we have 


to 4 decimal places. 


0.25 6 0.25 0.5 ee 
=, + 4e™-*° + 2e”-? + 4e®-’> + e?) 
0 


= 1.718 3188. 
(The true solution is. 1.718 281 8.) 
With h = 4 we have 


{ me oe + 4e%125 4 990.25 4 450.375 4 790.5 
3 | 
0 


i‘ 409-625 a. e975 she 4e9-875 rs e1) 
= 1.7182842 to 7 decimal places. 


This is an extremely accurate result. 


7. The approximate value for 


m/2 
L= { (1 + cos? x)!/? dx 
0 


cae re. 
using the composite Simpson’s method with h = : 1S 


Lx ‘c + cos?0)"/? + 411 + cos?— - 
aa 8 


1/2 3 1/2 
+ 21 + cos? + (1 + cos? 


sys 
+ f + cos? t 


= 1.910141 to 6 decimal places. 


(A more accurate answer, obtained by computer, using 
Simpson’s method with h = 7/200, is 1.910099 to 6 decimal 
places. The solution obtained is thus correct to 4 decimal 
places.) 


Solutions to the exercises in Section 5 


1. (i) From the hint we know that 


2 tan x sec? x = 2 tan x (1 + tan? x) 
= 2 tan x + 2tan* x. 


Since tan x is positive'and increasing in the interval 
0<x< 0.1, the maximum value of the above function occurs 
at x = 0.1. Hence, — 


2 tan c sec? c < 2 tan 0.1 sec” 0.1 
= 0.20 


(ii) The linear Taylor polynomial about 0 is given by 
p(x) =f (0) + xf'(0). 


For f(x) = tan x we have 


to 2 decimal places. 


(x) = sec? x = ——— 
f(x) as PA 


so that f(0) = 0 and f’(0) = 1. 


Hence the linear Taylor polynomial about 0 approximating 
tan x is given by 


p(x) = x. 
The error function is given by 
e(x) = —4x’f"(c,) 


where c,, lies between 0 and x. 
” d 2 rf 
Now f"(x) = 7 eed x 2 tan 906". x. 
x 


Hence 0 < f"(c,) = 2 tan c, sec? c, 
< 0.20 by part (i), since 0 < c, < 0.1. 
A bound for e(x) is therefore 
\e(x)| < 4x? x 0.20 = 0.1x?. 


0 1 
FAQS. 
1.222 1303 


1.1105 


2. The Newton-Raphson method uses the recurrence 
relation | 


2 ee 
ee 
We have 
f(x)=e& —3x, f'(x)=e& —3. 


6) 02 04°06 606 40-12 16> -z 
Figure 3. Graphs of y = e* and y = 3x. 


0 
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Hence for this problem the recurrence relation is 


exp (x,) — 3x, 


Xpt41 =X, — : 
es exp (x,) — 3 


(1) 


To obtain an initial approximation we sketch the graphs of 
y = e& and y = 3x as in Figure 3. 


From the graph we can see that two curves intersect at 
approximately x = 0.6 and x = 1.5. With x, = 0.6, Equation 
(1) above gives 


x, = 0.618 778 47 
x, = 0.619 061 22 
x3 = 0.619 061 29 
= X4q4. 
With xo = 1.5, Equation (1) gives 
x, = 1.5123361 
x, = 1.512 1346 


= X3. 


Hence to 5 decimal places the two roots are 0.61906 and 
A213. 


3. The Taylor series method of order 2 uses the recurrence 
relation 


¥41=¥,+hyY, + 5h’ Y,”, where 
y =*x7,4+1 and 
l= i > Xe, - 


With h = 0.1 we obtain the following table of results. 


Y= ¥%+%,¥, 


1 
1.21605 


Thus: when x = 0.1, y ~ 1.105, 
when x = 0.2, y ~ 1.222. 


4. (i) Let 
p(x) =a + bx 


Yui2=Y¥+01Y +0005), 


1.105 
1.222 1303 


be the interpolating polynomial. 


Since the graph of y = p(x) passes through the points (1, 1) 
and (24, 14), we have 


pil) =a+b=1 
p(24) = a + 24b = 14. 
Solving for a and b gives a = 3 and b = %. Hence 
p(x) = 0.6 + 0.4x. 
At x = 1.44 we have 
p(1.44) = 1.176. 


(ii) Since }x~%/? is a positive decreasing function on the 

interval 1 < x < 24, its maximum value is at x = 1, and is 

equal to 4. Hence 
be {<4 fori<cem% 

(iii) From Theorem 4.1 on page 39, the error function is 

given by 


e(x) = —(x — 1)(x — 24) 


Now f (x) = mes and so 
f'(x) =4x7'2 — and 


I"(cx) 
= 


f"(x) = —x 9? 
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From part (11) we know that an upper bound for | f”(c,.)| 
where 1 < c, < 23 is given by M = 4. Hence at x = 1.44, 


|(1.44 — 1)(1.44 — 2.25) 
8 


The actual error at x = 1.44 is 
£(1.44) = p(1.44) — ./1.44 = 1.176 — 1.2 
— 0.024 


and its magnitude is less than the error bound for the 
approximation. 


|e(1.44)| < = 0.044 55. 


5. Simpson’s method in composite form with h = 0.25, 


approximating the integral of tan x between 0 and 1, gives 


: 0.25 
{ tan x dx ~ —z (tan 0+ 4 tan 0.25 + 2 tan 0.5 
0 


+ 4tan 0.75 + tan 1) 
= 0.616 480 52. 


55 


(The true solution is 
1 
{ tan x dx = [log, (sec x)]§ = 0.615 626 47.) 
0 
6. From Section 3 the length L of a curve y is given by 
d 2 
L=[ hi + fa dx. 
me dx 
dy 
For y =e" we have —_ = e* and 


1 
L=| 1 +e" dx. 


0 
Simpson’s method for this integral, with h = 0.25, gives 


b= /2 +41 +72 42/14e4+4/1 + 23? 
+ ./1 + e?) = 2.003 5275. 


7 
20 nN eke oo Oo nbs nad enses eee 
2000 RC8888 6) 2oeeee ea seee sae 
20g 088el ooo seems oo enaeas 
—S88G0E 98800) 665 ent oeenee 
ao ou 
“2 eeneOsU esos e2ee a" 
leet eeeeeeC re Bnanene 
“6 00Ce OneU eee sUeoosseene 
eT 
Seen GGGOEe hen R ee eoeee ie 
“0 0G sGtth see eeeeee eee 
| OSG nee es eee eee eee aG 
So CeRORe ee es Cees ee eesaGEe 
2G U008 ee GeeOn neon aees 
J teeCe® onl eetneoeese 
20 ee eee eeeeOeeans ene senes 
ones [ee0nee0 2oeceE 
me eOGke eG te eens eee 
"ECS shen s9s enero seen eae ess 
C00 208 26RKr Chee Annee Ceenens 
Tye eee peer 
Ti a oid Gado 6 
TS EEE cceuueue4o See eosens 
ogee eee eae 


om ou 


